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Abstract. We present the third in the series of papers describing Poisson properties of pla- 
nar directed networks in the disk or in the annulus. In this paper we concentrate on special 
networks N„ ,, in the disk that correspond to the choice of a pair (u, v) of Coxeter elements 
in the symmetric group S„ and the corresponding networks N°y in the annulus. Boundary 
measurements for N,, ,. represent elements of the Coxeter double Bruhat cell G" ' c GL„. 
The Cartan subgroup H acts on G"''' by conjugation. The standard Poisson structure on 
the space of weights of N„_,, induces a Poisson structure on G" ' , and hence on the quo- 
tient G"'*"/!!, which makes the latter into the phase space for an appropriate Coxeter-Toda 
lattice. The boundary measurement for N° ,, is a rational function that coincides up to a 
nonzero factor with the Weyl function for the boundary measurement for Nu,v The corre- 
sponding Poisson bracket on the space of weights of N° induces a Poisson bracket on the 
certain space K,, of rational functions, which appeared previously in the context of Toda 
flows. 

Following the ideas developed in our previous papers, we introduce a cluster algebra 
^ on K„ compatible with the obtained Poisson bracket. Generalized Backlund-Darboux 
transformations map solutions of one Coxeter-Toda lattice to solutions of another preserv- 
ing the corresponding Weyl function. Using network representation, we construct gener- 
alized Backlund-Darboux transformations as appropriate sequences of cluster transforma- 
tions in y(. 



1. Introduction 

This is the third in the series of papers in which we investigate Poisson geometry of 
directed networks. In [21, 22], we studied Poisson structures associated with weighted 
directed networks in a disk and in an annulus. The study was motivated in part by Poisson 
properties of cluster algebras. In fact, it was shown in [21] that if a universal Poisson 
bracket on the space of edge weights of a directed network in a disk satisfy an analogue of 
the Poisson-Lie property with respect to concatenation, then the Poisson structure induced 
by this bracket on the corresponding Grassmannian is compatible with the cluster algebra 
structure in the homogeneous coordinate ring of the Grassmannian. In this paper we deal 
with an example that ties together objects and concepts from the theory of cluster algebras 
and directed networks with the theory of integrable systems. 

Integrable systems in question are the Toda flows on GL„. These are commuting Hamil- 
tonian flows generated by conjugation-invariant functions on GL„ with respect to the stan- 
dard Poisson-Lie structure. Toda flows (also known as characteristic Hamiltonian systems 
[30]) are defined for an arbitrary standard semi-simple Poisson-Lie group, but we will 
concentrate on the GL„ case, where as a maximal algebraically independent family of 
conjugation-invariant functions one can choose Fk '■ GL„ 3 X jtrX'', k = l,...,n-l. 
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The equation 
(1.1) 



of motion generated by Fk has a Lax form: 

X. -i(»4X*)-/r-(X*)) 



dt 



where n+{A) and n^iA) denote strictly upper and lower parts of a matrix A. 

Any double Bruhat cell G"'^, m, v g 5 „, is a regular Poisson submanifold in GL„ invariant 
under the right and left multiplication by elements of the maximal torus (the subgroup of 
diagonal matrices) H c GL„. In particular, G"'^ is invariant under the conjugation by 
elements of H. The standard Poisson-Lie structure is also invariant under the conjugation 
action of H on GL„. This means that Toda flows defined by (1.1) induce commuting 
Hamiltonian flows on G"'^7H where H acts on G"'^ by conjugation. In the case when 
V = M ' - {n 12 . . .n - V), G"'^ consists of tridiagonal matrices with nonzero oflf-diagonal 
entries, C'^'/H can be conveniently described as the set Jac of Jacobi matrices of the form 

foi 1 ••• 
ai bi 1 ••• 



(1.2) 



L = 







1 

bn 



ai---a„-i 7^0, detL^tO. 



Lax equations (1.1) then become the equations of the finite nonperiodic Toda hierarchy: 



dt 



the first of which, corresponding to A: = 1, is the celebrated Toda lattice 
d 



dt 

dt^^ 



= ajibj+i - bj), 
= (aj-aj-i), 



j=l,...,n-l, 



with the boundary conditions ao = a„ = 0. Recall that detL is a Casimir function for the 

standard Poisson-Lie bracket. The level sets of the function det L foliate Jac into 2(« - 1)- 
dimensional symplectic manifolds, and the Toda hierarchy defines a completely integrable 
system on every symplectic leaf. Note that although Toda flows on an arbitrary double 
Bruhat cell G"'*' can be exactly solved via the so-called factorization method (see, e.g. 
[31]), in most cases the dimension of symplectic leaves in G"'" jM. exceeds 2(n - 1), which 
means that conjugation-invariant functions do not form a Poisson commuting family rich 
enough to ensure Liouville complete integrability. 

An important role in the study of Toda flows is played by the Weyl function 



(1.3) 



w(A) = m{A;X) = ((Al - X)-'ei, d) 



where p(A) is the characteristic polynomial of X and q(A) is the characteristic polynomial 
of the (n - 1) X (n - 1) submatrix of X formed by deleting the first row and column (see, e.g., 
[7, 28, 5]). Difl'erential equations that describe the evolution of miA;X) induced by Toda 
flows do not depend on the initial value X(0) and are easy to solve: though nonlinear, they 
are also induced by linear differential equations with constant coefficients on the space 



(1.4) 



\M{A) = 



by the map M{X) i-> m{A) 



: degP = n, degQ = n - 1, P, 2 are coprime, P(0) Oj 
: Um^_c« AM{A) + 0. 



-^M(-/l), where //o 
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It is easy to see that m(A; X) is invariant under the action of H on G"'^' by conjugation. 
Thus we have a map from C'^/H into the space 



In the tridiagonal case, this map, sometimes called the Moser map, is invertible: it is a 
classical result in the theory of moment problems that matrix entries of an element in 
Jac can be restored from its Weyl function m{A; X) via determinantal formulas for matrix 
entries of X in terms of Hankel determinants built from the coefficients of the Laurent 
expansion of m(/l; X). These formulas go back to the work of Stieltjes on continuous 
fractions [33] (see, e.g. [1] for details). 

In this paper, we study double Bruhat cells G"'^ that share common features with the 
tridiagonal case: 

(i) the Toda hierarchy defines a completely integrable system on level sets of the deter- 
minant in G"'''/H, and 

(ii) the Moser map m„,v : G"'''/H — > "Wn defined in the same way as in the tridiagonal 
case is invertible. 

We will see that double Bruhat cells G"'" associated with any pair of Coxeter elements m, v e 
S n enjoy these properties. (Recall that a Coxeter element in 5 „ is a product of n - 1 distinct 
elementary transpositions.) Double Bruhat cells of this kind has previously appeared (for 
an arbitrary simple Lie group) in [23] in the context of integrable systems and in [3, 36] in 
connection with cluster algebras of finite type. We will call any such double Bruhat cell 
a Coxeter double Bruhat cell. Integrable equation induced on G"'^/H by Toda flows will 
be called Coxeter-Toda lattices. This term was first used in [23] in the case u = v for an 
arbitrary simple Lie group, which generalizes the relativistic Toda lattice that corresponds 
to the choice u = v= s„_i • • • si in GL„. In [12, 13], the corresponding integrable systems 
for V - Sn-\ ■ ■ ■ si and an arbitrary Coxeter element u were called elementary Toda lattices. 
In the latter case, G"'"/!! can be described as a subset of Hessenberg matrices subject to 
certain rank conditions on submatrices. The tridiagonal case corresponds to the choice 

V = Sn-\ ■ ■ ■ SuU - Sx - ■ ■ Sn-\. 

Since Coxeter-Toda flows associated with diflferent choices of (m, v) lead to the same 
evolution of the Weyl function, and the corresponding Moser maps are invertible, one 

can construct transformations between different G"'^ /H that preserve the corresponding 
Coxeter-Toda flows and thus serve as generalized Bdcklund-Darboux transformations be- 
tween them. 

Our goal is to describe these transformations from the cluster algebra point of view. To 
this end, we construct a cluster algebra of rank 2n-2 associated with an extension of the 
space (1.4) 



(Note that ^„ is embedded into !R„ as a codimension 1 subspace.) Distinguished clusters 
x„ y in this algebra correspond to Coxeter double Bruhat cells, and are formed by certain 
collections of Hankel determinants built out of coefficients of the Laurent expansion of 
an element in Sequences of cluster transformations cormecting these distinguished 
clusters are then used as the main ingredient in the construction of generalized Backlund- 
Darboux transformations. 

The insight necessary to implement this construction is drawn from two sources: 



^V„ = \m(A) = 



q(A) 



: degp = n, degq = n- 1, p, ^ are monic and coprime, p(0) ^ 
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(i) the procedure for the inversion of the Moser map, that can be viewed as a general- 
ization of the inverse moment problem, and 

(ii) interpretation of functions in as boundary measurement functions associated with 
a particular kind of networks in an annulus. 

Before discussing the organization of the paper, we would like to make two remarks. 
First, birational transformations between G"'""-^'"^^ /H and G"''""-^'"^^ /H for two different 
Coxeter elements u, u' that serve as generaUzed Backlund-Darboux transformation be- 
tween the corresponding elementary Toda lattices were first studied in [12]. Second, a 
cluster algebra closely related to the one we considered here recently appeared in [25] and 
was a subject of a detailed combinatorial study in paper [9], where cluster mutations along 
the edges of a certain subgraph of its exchange graph were shown to describe an evolution 
of an An type Q-system — a discrete evolution that anses in the analysis of the XXX-model, 
which is an example of a quantum integrable model. In [9], solutions of the Q-system are 
represented as Hankel determinants built from coefficients of a certain generating function, 
that turns out to be rational and can be represented as a matrix element of a resolvent of an 
appropriate linear operator. 

The paper is organized as follows. 

In Section 2 we go over the necessary background information on double Bruhat cells, 
Toda flows, cluster algebras, networks and associated Poisson structures. We then proceed, 
in Section 3, to describe a parametrization of a Coxeter double Bruhat cell. This is a 
particular case of the Berenstein-Fomin-Zelevinsky parametrization [2, 15]: for a generic 
element X in G"'^, we consider a factorization of X into elementary bidiagonal factors 
consistent with the Gauss factorization of X, that kX - X^XqX+, where Xq is the diagonal 
matrix diag(t/i, . . . , t/„), X+ is the product of n - 1 elementary upper bidiagonal factors 
Et{ct), ! = 1 , . . . , n - 1 , with the order of factors in the product prescribed by v, and X- is 
the product of n - 1 elementary lower bidiagonal factors Ej{cJ), i = \, . . . ,n - \, with the 
order of factors in the product prescribed by u. We also give an intrinsic characterization 
of a double Bruhat cell. 

Elements G"'"/!! are parametrized by d, and c,- = c^c~, j = 1, . . . ,n - 1. In Section 
4 we show that these parameters can be restored as monomial expressions in terms of an 
appropriately chosen collection of Hankel determinants built from the coefficients of the 
Laurent expansion of the Weyl function m{X). (In [14], a similar inverse problem was 
solved for the case v = Sn-i •■• si,u arbitrary.) Both the choice of Hankel determinants and 
exponents entering monomial expressions for di, ci are uniquely determined by the pair 
(m, v). 

In Section 5, the map X i-> m{A\X) is given a combinatorial interpretation in terms of 
weighted directed planar networks. To an elementary bidiagonal factorization of X e G"'" 
there corresponds a network Nu,v in a square (disk) with n sources located on one side of 
the square and n sinks located at the opposite side, both numbered bottom to top (see, e.g. 
[15, 16, 1 1]). By gluing opposite sides of the square containing sinks and sources in such 
a way that each sink is glued to the corresponding source and adding two additional edges, 
one incoming and one outgoing, one obtains a weighted directed network in an annulus 
(outer and inner boundary circles of the annulus are formed by the remaining two sides of 
the square). Networks in an annulus were studied in [22]. The network we just described, 
A^°y, has one sink and one source on the outer boundary of an annulus and, according to 
[22], the boundary measurement that corresponds to this network is a rational function 
M{A) in an auxiliary parameter A. We show that -M{-A) is equal to m{A\X) times the 
product of weights of the incoming and outgoing edges in A'"^. 
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The determinantal formulae for the inverse of the Moser map are homogeneous of de- 
gree zero with respect to coefficients of the Laurent expansion, therefore the same formulae 
applied to M(A) also recover c„ d,. Thus, we can define a map p„_v : (C*)^" -> G"'^/K in 
such a way that the through map 

is the identity map. 

In the remainder of Section 5, we use the combinatorial data determined by the pair 
(m, v) (or, in a more transparent way, by the corresponding network A^°^) to construct a 
cluster algebra J?l = ~9lu,v, with the (slightly modified) collection x„_y serving as the ini- 
tial cluster. The matrix fi„_v that determines cluster transformations for the initial cluster 
is closely related to the incidence matrix of the graph dual to N°^^. To construct ^u,v, we 
start with the Poisson structure induced on boundary measurement functions by a so-called 
standard Poisson bracket on the space of face weights associated with ^ (this bracket 
is a particular case of the general construction for networks in the annulus given in [22J). 
Initial cluster variables, viewed as functions on form a coordinate system in which 
this Poisson structure takes a particular simple form: the Poisson bracket of logarithms 
of any two functions in the family is constant. This allows us to follow the strategy from 
[20] to construct ^u,v as a cluster algebra compatible with this Poisson bracket. We then 
show that ^u,v does not depend on the choice of Coxeter elements u, v, that is, that for any 
(«', v'), the initial seed of ^u'y is a seed in the cluster algebra ^u,v Therefore, the change 
of coordinates T^ f : x„_v Xu'y is accomplished by a sequence of cluster transforma- 
tions. Moreover, the ring of regular functions on Kn coincides with the locaUzation of the 
complex form of J?l with respect to the stable variables. We complete Section 5 with the 
discussion of the interplay between our results and those of [9]. In particular, we provide 
an alternative proof for one of the main results of [9] concerning the Laurent positivity of 
the solutions of Q-systems. 

In the final section, we interpret generalized Backlund-Darboux transformations be- 
tween Coxeter-Toda lattices corresponding to diflTerent pairs of Coxeter elements in terms 
of the cluster algebra J?l by observing that the map 

(1.5) = p,,y o r«>' o T„,, : G«-7H ^ G^'-^'/H, 

with T„_v being the right inverse of p„_v. preserves flows generated by conjugation-invariant 
functions and makes the diagram 

G"'7H ^ G« /H 




commutative. We obtain explicit formulas for tr" ;.' and, as a nice application, present for- 
mulas that transform solution of the usual Toda lattice into solutions of the relativistic one. 
Besides, we explain how one represents generaUzed Backlund-Darboux transformations 
as equivalent transformations of the network N° ,,. Finally we show that classical Darboux 
transformations are also related to cluster algebra transformations via a formula similar 
to (1.5). 
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2. Preliminaries 

In this section we collect the necessary background information on double Bruhat cells, 
Toda flows and directed networks on surfaces. Though notions and results that we will 
need on the first two subjects can be as easily stated for an arbitrary semisimple group, we 
will limit ourselves to the GL„ case. 

2. 1 . Double Bruhat cells. Let b+, n+, b_, n_ be, resp., algebras of upper triangular, strictly 
upper triangular, lower triangular and strictly lower triangular matrices. 

The connected subgroups that correspond to b+, b_, Tt+, Tt_ will be denoted by B+, B_, 
N+, N_. We denote by H the maximal torus (the subgroup of diagonal matrices) in GL„. 

Every ^ e gl„ can be uniquely decomposed into 

^ = ^_ + fo + ^+, 

where f+ e n+, ^_ E n_ and is diagonal. Consequently, for every X in an open Zariski 
dense subset of GL„ there exists a unique Gauss factorization 

X = X.XoX+, Z+ e N+, Z_ e N_, Xq e H. 

Let s,-, / 6 [1, « - 1], denote the elementary transposition (;, ; + 1) in the symmetric group 
5„. A reduced decomposition of an element w G 5„ is a representation of w as a product 
w = i„ • • • of the smallest possible length. A reduced decomposition is not unique, but 
the number I depends only on w and is called the length of w and denoted by Kyv). The 
sequence of indices i = (i'l, . . . , ii) that corresponds to a given reduced decomposition of w 
is called a reduced word for w. The unique element of 5 „ of maximal length (also called 
the longest element of S „) is denoted by wq. 

We will also need need a notion of a reduced word for an ordered pair (m, v) of elements 
in ^n. It is defined as follows: if (j'l, . . . , j;(„)) is a reduced word for u and {i\, i'^^^-^ is 
a reduced word for v, then any shuflle of the sequences {i\,. . ., ii(u)) and . . . , -i'i(y^) is 
called a reduced word for (w, v). 

Let us fix an embedding ofS„ into GL„ and denote the representative of w € 5 „ in GL„ 
by the same letter w. The Bruhat decompositions of GL„ with respect to B+ and B_ are 
defined, resp., by 

GL„ = U„es„B+MB+, GL„ = Uves„B_vB_. 
For any m, v e the double Bruhat cell is defined as 

G''^ = B+MB+nB_vB_. 

According to [15], the variety G"'^ is biregularly isomorphic to a Zariski open subset of 

i^;(«)+;(v)+n ^ corresponding birational map from c'(")+'(i')+« j-q qu.v ^.^^ constructed quite 

explicitly, though not in a unique way. Namely, fix a reduced word i for the pair (u, v) and 
consider, in addition, a sequence k = {k\,. ..,k„) obtained as an arbitrary re-arrangement 
of numbers t through in. Let j = (j'l , . . . , ji(u)+i(y)+n) be a shuffle of k and i; we set Oiji) to 
+ if ji > 0, to - if ji < 0, and to if ji e k. 

Denote by etj an elementary nx n matrix (<5ia<J;)s)^^_j- For t € C, i,j € [l,n - I], 
A;e[l,n], let 

(2.1) E;it) = l + teM,i, E^it) = l + tejj^u E°it) = 1 + (t - l)ekk. 
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Then the map Xj : c'(">+'(^)+'' ^ G"'" can be defined by 

;(«)+;(v)+n 

(2.2) Xj(t)= f] E^^'a,). 

9=1 

Parameters , . . . , tnu)+i(v)+n constituting t are called factorization parameters. ExpUcit 
formulae for the inverse of the map (2.2) in terms of the so-called twisted generalized 
minors were found in [15]. 

2.2. Toda flows. Next, we review the basic facts about the Toda flows on GL„. 
Recall that the standard Poisson-Lie structure on GL„ is given by 

{/i,/2}5L„(X) = ^(/?(V/i(X)X),V/2(X)Z)- i(/?(XV/i(X)),XV/2(X)), 

where ( , ) denotes the trace-form, V is the gradient defined with respect to the trace-form, 
and R : gl„ gl„ is the standard R-matrix given by 



/?(^=^,-^_ = (sign(7-0^,v)".^j 



Double Bruhat cells are regular Poisson submanifolds of GL„ equipped with the stan- 
dard Poisson-Lie structure (see [30, 26, 35]). Furthermore, 

(i) any symplectic leaf of GL„ is of the form 5 "'"a, where S"'^ c G"''' is a certain 
distinguished symplectic leaf and a is an element of the Cartan subgroup, and 

(ii) the dimension of symplectic leaves in G"-^ equals /(m) -i- /(v) + corank(Mv"' - Id), see 
[30, 26]. 

Conjugation-invariant functions on GL„ form a Poisson-commuting family (see, e.g., 
[31]). Any such function F generates a Hamiltonian flow described by the Lax equation 



(2.3) dX/dt = 



X, -i^(XVF(X)) 



The resulting family of equations is called the hierarchy of Toda flows (in [30], the term 
characteristic Hamiltonian systems is used). If one chooses F(X) = Fk(X) = i tr(X*), 
then equation (2.3) becomes (1.1). Functions F\,..., F„_i form a maximal family of alge- 
braically independent conjugation-invariant functions on GL„. 

For an element h e GL„, denote by Ch the action of h on GL„ by conjugation: ChiX) = 
hXh~^. For any smooth function / on GL„ we have 

V(/oC,)=A4-i(V/). 

Furthermore, if h belongs to H, then it is easy to see that 

RiAdhA^))^Adh->(R{^)) 

for any ^ e gl„. Together, these observations imply that for any /i e H and any pair of 
smooth functions /i, /i on GL„, 

{/l°Q,/2oC,} = {/i,/2}oCft. 

In other words, the action of H on GL„ by conjugation is Poisson with respect to the stan- 
dard Poisson-Lie structure. Since the action preserves double Bruhat cells, the standard 
Poisson-Lie structure induces a Poisson structure on G"'''/H, and the Toda hierarchy in- 
duces the family of connmuting Hamiltonian flows on G"'^/H. 
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Remark 2.1. (i) The Lax equation (2.3) can be solved explicitly via the factorization 
method [3 1], which we will not review here. 

(ii) Written in terms of matrix entries, equations (2.3) have exactly the same form as 
equations of the Toda hierarchy on gin, where the relevant Poisson structure is the Lie- 
Poisson structure associated with the R-matrix Lie bracket 7/]^? = ^ ([/?(^), 77] + /?(??)]). 
In fact, viewed as equations on the algebra of n x n matrices, the Toda hierarchy becomes 
a family of bi-Hamiltonian flows with compatible hnear and quadratic Poisson brackets 
given by, respectively, Lie-Poisson and the extension of the Poisson-Lie brackets. How- 
ever, we will not need the linear Poisson structure in the current paper. 

2.3. Cluster algebras and compatible Poisson brackets. First, we recall the basics of 
cluster algebras of geometric type. The definition that we present below is not the most 
generid one, see, e.g., [17, 3| for a detailed exposition. 

The coefficient group P is a free multiplicative abelian group of a finite rank m with 
generators gi,. . . ,gm- An ambient field is the field g of rational functions in n indepen- 
dent variables with coefficients in the field of fractions of the integer group ring ZP = 
1-{g\^, . . . , g^} (here we write x^^ instead of x, x"'). It is convenient to think of 5 as of the 
field of rational functions in n + m independent variables with rational coefficients. 

A seed (of geometric type) in g is a pair 2 - (x, B), where x = (xi , . . . , Xn+m), xi, . . . ,x„ 
is a transcendence basis of g over the field of fractions of ZP, x„+, - gj for / e [1, m], and 
5 is an n X (n -H m) integer matrix whose principal part (that is, the « x n submatrix formed 
by the columns 1, . . . , n) is skew- symmetric. The (« + OT)-tuple x is called a cluster, its 
elements called cluster variables, and x„+\, x„+m are stable variables. 

Given a seed as above, the cluster transformation in direction k€ [1, n] is defined by 

X Xt = (X \ [Xk]) U {Xk], 

where the new cluster variable Xk is given by the exchange relation 
(2.4) XkXk^ Y] n 

\<i<n+m \<i<n+m 
bn>0 bn<0 

here, as usual, the product over the empty set is assumed to be equal to 1. 
We say that B is obtained from B by a matrix mutation in direction k if 

{-bij, ifi = k or j = k; 

\bik\bkj + bik\bkj\ 
bij-\ , otherwise. 

Given a seed S = (x, B), we say that a seed S = (x, B) is adjacent to E (in direction k) 
if X is obtained from x and B is obtained from B by a cluster transformation and a matrix 
mutation, respectively, in direction k. Two seeds are mutation equivalent if they can be 
connected by a sequence of pairwise adjacent seeds. The cluster algebra (of geometric 
type) - !MB) associated with Z is the ZP-subalgebra of g generated by all cluster 
variables in all seeds mutation equivalent to 2. The complex form of J?l is defined as Ji 
tensored by C and is denoted J?lc. 

Let y be a Zariski open subset in C"^™, J?lc be the complex form of a cluster algebra 
of geometric type. We assume that the variables in some extended cluster are identified 
with a set of algebraically independent rational functions on V. This allows us to identify 
cluster variables in any cluster with rational functions on V as well, and thus to consider 
J?lc as a subalgebra of the field C(y) of rational functions on V. Finally, we denote by 
the locaUzation of J?lc with respect to the stable variables that do not vanish on V. 
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Proposition 2.1. Let V and J?l as above satisfy the following conditions: 

(i) each regular function on V belongs to J?!^; 

(ii) there exists a cluster x = {xi,..., x„+m) in Jic consisting of algebraically indepen- 
dent functions regular on V; 

(iii) any cluster variable x^, k e [l,n], obtained by the cluster transformation (2.4) 
applied to x is regular on V. 

Then is isomorphic to the ring 0(V) of regular functions on V. 

Proof. All we have to prove is that any element in is a regular function on V. The 
proof follows the proof of a similar statement for double Bruhat cells in [37] and consists 
of three steps. 

Lemma 2.1. Let z = (zi , . . . , z„+m) be an arbitrary cluster in J?Ic. If a Laurent monomial 

M — zf^ • ■ • Z^'+m regular on V then di > O for ; e [1, nj. 

Proof. Indeed, assume that <0 for some k€ [1 , n] and consider the cluster z^.. By (2.4), 
M can be rewritten as M = Miz^*//""*, where Mi is a Laurent monomial in common 
variables of z and z^, and P is the binomial (in the same variables) that appears in the right 
hand side of (2.4). By condition (i) and the Laurent phenomenon (Theorem 3.1 in [17]), 
M can be written as a Laurent polynomial in the variables of z^. Equating two expressions 
for M, we see that P"* times a polynomial in variables of z^ equals a Laurent monomial in 
the same variables. This contradicts the algebraic independence of variables in z^, which 
follows from the algebraic independence of variables in x. □ 

Lemma 2.2. Let zbe a cluster variable in an arbitrary cluster x, and assume that z is a 
regular function on V. Then z is irreducible in the ring of regular functions on V. 

Proof. Without loss of generality, assume that z„+\ = -V'n+i , . . . , z„+m' - Xn+m' do not vanish 
on V, and Zn+m'+i = x„+m'+i, . . .,Zn+m = Xn+m may vanish on V. Moreover, assume to the 
contrary that z = fg, where / and g are non-invertible regular functions on V. By condi- 
tion (i) and Proposition 11. 2 of [18], both / and g are Laurent polynomials in zi , . . . , z„+m' 
whose coefficients are polynomials in Zn+m'+i> ■ ■ ■,Zn+m- Applying the same argument as 
in the proof of Lemma 2.1, we see that both / and g are, in fact, Laurent monomials in 
Zi , . . . , Zn+m and that Zn+i , . . . , z«+m' enter both / and g with a non-negative degree. More- 
over, by Lemma 2.1, each cluster variable zi , . . . , z„ enters both / and g with a non-negative 
degree. This can only happen if one of / and g is invertible in 0{V), a contradiction. □ 

Denote by Uq aV the locus of aU ? € V such that x,(?) 9^ for aU j e [1, «]. Besides, 
denote hyUk<zV the locus of aU ? e V such that Xi{t) for aU j € [1, n] \ A: and Xk{t) + 0. 

Lemma 2.3. Let U = U^^gf/,, then codimV \ U >2. 

Proof Follows immediately from Lemma 2.2 and conditions (ii) and (iii). □ 

Assume that there exists / e that is not regular on V. Recall that V is the comple- 
ment of a finite union of irreducible hypersurfaces D, in C"^™. Therefore, the divisor of 
poles of / has codimension 1 in C"^™. Since / is not regular on V, this latter divisor does 
not lie entirely in the union of Di, and hence its intersection with V has codimension 1 in 
V. Therefore, by Lemma 2.3, it intersects U nontrivially. To complete the proof, note that 
by Proposition 1 1.2 of [18], any function in is regular on U, a contradiction. □ 

Let {•, •} be a Poisson bracket on the ambient field 5. We say that it is compatible with 
the cluster algebra J?l(B) if, for any cluster x = {xi,..., x„+m), one has (x„ xj] = (OijXiXj, 
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where ojij e Z are constants for all i,] e [\,n + m\. The matrix = is called 
the coefficient matrix of {-, •) (in the basis x); clearly, is skew- symmetric. A complete 
description of Poisson brackets compatible with in the case rankB = n is given in 
[20]. 

2.4. Networks on surfaces with boundaries. Let 5 be a disk with c > holes, so that 
its boundary dS has c + 1 connected components, and let G = (V, E) be a directed graph 
embedded in S with the vertex set V and the edge set E. Exactly r of its vertices are 
located on the boundary dS . They are denoted bi,...,br and called boundary vertices. 
Each boundary vertex is labeled as a source or a sink. A source is a vertex with exactly 
one outcoming edge and no incoming edges. Sinks are defined in the same way, with 
the direction of the single edge reversed. The number of sources is denoted by n and the 
number of sinks hy m = r - n. All the internal vertices of G have degree 3 and are of two 
types: either they have exactly one incoming edge, or exactly one outcoming edge. The 
vertices of the first type are called (and shown on figures) white, those of the second type, 
black. 

A pair (v, e), v e V, e e E, is called a flag if v is an endpoint of e. To each flag (v, e) we 
assign an independent variable Xy.^. Let u and v be two endpoints of e. The edge weight 
We is defined by Wg = Xv,eXu,e- A perfect network N = (G, w,pi, . . . ,pe) is obtained from 
G weighted as above by adding c nonintersecting oriented curves p, (called cuts) in such 
a way that cutting S along p, makes it into a disk (note that the endpoints of each cut 
belong to distinct connected components of dS). The points of the space of edge weights 
£>N = (R\0)'^' (or (C\0)'^') can be considered as copies of the graph G with edges weighted 
by nonzero numbers obtained by speciahzing the variables Xv,e to nonzero values. 

Assign an independent variable /I, to each cut p;. The weight of a path P between two 
boundary vertices is defined as the product of the weights of all edges constituting the 
path times a Laurent monomial in /I,. Each intersection point of P with p, contributes to 
this monomial /l, if the oriented tangents to P and p, at this point form a positively oriented 
basis, and Aj^ otherwise (assuming that all intersection points are transversal). Besides, the 
sign of the monomial is defined via the rotation number of a certain closed curve built from 
P itself, cuts and arcs of dS . For a detailed description of the corresponding constructions, 
see [21, 29] in the case c - (networks in a disk, no cuts needed, the path weight is a 
signed product of the edge weights) and [22] in the case c - I (networks in an annulus, 
one cut p and one additional independent variable A involved, the path weight is a signed 
product of the edge weights times an integer power of A). The boundary measurement 
between a source bi and a sink bj is then defined as the sum of path weights over all (not 
necessary simple) paths from bi to bj. It is proved in the above cited papers that a boundary 
measurement is a rational function in the weights of edges (in case of the disk) or in the 
weights of edges and A (in case of the annulus). 

Boundary measurements are organized in the boundary measurement matrix, thus giv- 
ing rise to the boundary measurement map from to the space of n X m matrices (for 
c = 0), or the space of nxm rational matrix functions (for c = 1). The gauge group acts on 
Sn as follows: for any internal vertex v of A'^ and any Laurent monomial L in the weights 
We of N, the weights of all edges leaving v are multiplied by L, and the weights of all edges 
entering v are multiplied by L Clearly, the weights of paths between boundary vertices, 
and hence boundary measurements, are preserved under this action. Therefore, the bound- 
ary measurement map can be factorized through the space 'Fn defined as the quotient of 
&N by the action of the gauge group. In [22] we explained that T^n can be identified with 
the relative cohomology group H^{G, GC\dS) with coefficients in the multipUcative group 
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of nonzero real numbers. This gives rise to the representation 

Tn = //'(G U dS)IHHdS) e H^\dS)IH\G UdS)^r^(Br^ 



The space T^j^ can be described as follows. The graph G divides S into a finite number 
of connected components called /ace*. The boundary of each face consists of edges of G 
and, possibly, of several arcs of 55 . A face is called hounded if its boundary contains only 
edges of G and unbounded otherwise. Given a face /, we define its face weight y/ as the 
function on £>i^ that assigns to the edge weights We,e&E, the value 



where = 1 if the direction of e is compatible with the counterclockwise orientation of 
the boundary df and 7e = ~1 otherwise. Face weights are invariant under the gauge group 
action, and hence are functions on Tj^, and, moreover, form a basis in the space of such 
functions. 

In [21, 22] we have studied the ways to turn the space of edge weights into a Poisson 
manifold by considering Poisson brackets on the space of flag variables satisfying certain 
natural conditions. We proved that all such Poisson brackets on form a 6-parameter 
family, and that this family gives rise to a 2-parameter family of Poisson brackets on fj^. 
In what follows we are interested in a specific member of the latter family (obtained by 
setting or = 1/2 and ;S = -1/2 in the notation of [21, 22]). For the reasons that wiU be ex- 
plained later, we call this bracket the standard Poisson bracket on f^- The corresponding 
4-parameter family of Poisson brackets on &n is called standard as well. 

Given a perfect network A' as above, define the directed dual network N* = (G*, w*) 
as follows. Vertices of G* are the faces of A'^. Edges of G* correspond to the edges of 

that connect either two internal vertices of different colors, or an internal vertex with a 
boundary vertex; note that there might be several edges between the same pair of vertices 
in G* . An edge e* of G* corresponding to e is directed in such a way that the white endpoint 
of e (if it exists) lies to the left of e* and the black endpoint of e (if it exists) lies to the right 
of e. The weight w*(e*) equals 1 if both endpoints of e are internal vertices, and 1/2 if one 
of the endpoints of e is a boundary vertex. 

Proposition 2.2. The restriction of the standard Poisson bracket on Ttf to the space Tj^ is 
given by 



For networks in a disk, the above proposition is a special case of Lemma 5.3 of [21]. 
For other surfaces the proof is literally the same. 

In what follows, we will deal with networks of two kinds: acyclic networks in a disk 
with the same number of nonalternating sources and sinks, and networks in an annulus 
obtained from the networks of the first kind by a certain construction, to be described 
below. 

In the former case we assume that n sources are numbered clockwise and are followed 
by n sinks numbered counterclockwise. The weight of a path in this case is exactly the 
product of edge weights involved. The boundary measurements are organized into a n x 
n matrix X in such a way that Xij is the boundary measurement between the ;th source 
and jth sink. One can concatenate two networks of this kind by gluing the sinks of the 
former to the sources of the latter. If Xi, X2 are the matrices associated with the two 
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networks, then the matrix associated with their concatenation is XiX2. This fact can be used 
to visualize parametrization (2.2). Indeed, an n x n diagonal matrix diag(c/i, ...,dn) and 
elementary bidiagonal matrices Et(1) and Ej{u) defined by (2.1) correspond to building 
blocks shown on Figure 1 a, b and c, respectively; all weights not shown explicitly are 
equiil to 1. Note that building blocks themselves are not networks, since their edge weights 
do not comply with the rules introduced above. However, as we will see below, objects 
glued from building blocks comply with all the rules. 



' i+l j+l • ^ — j+l 

2 . ^ 2 i . ^ ^ . j .^J—^ j 

i 

a) b) c) 



Figure 1 . Three building blocks used in matrix factorization 

The concatenation of n(n - 1) building blocks of the second and the third types and one 
building block of the first type, in an appropriately chosen order and with each building 
block having its own nontrivial weights, describes a generic element of GL„ (see, e.g. [11]). 
The structure of the obtained network is given by Figure 2. Here and in what follows we 
use the gauge group action to decrease the number of parameters of networks in question. 
In particular, this network has 2n(n - 1) internal vertices, therefore, one can use the gauge 
group action to change the weights of 2n(n - 1) edges to 1. It is convenient to choose 
these edges to be all the horizontal edges except for one middle edge in each horizontal 
chain. The weights on the remaining edges are Laurent monomials in the initial weights of 
the network. For example, if the endpoints of an edge e belong to levels / and i+l, then 
Ue = WeWp.^^ /wp^, where f, and Pi+i are the horizontal paths from the endpoints of e to the 
sinks i and i+l, respectively. 




Figure 2. Generic planar network; the weights of edges are Laurent 
monomials of the initial weights 

The following result, which is a special case of Theorem 4.1 from [21], explains why 
we call the bracket in consideration standard. 

Theorem 2.1. For any network N as above with n sources and n sinks the map from 
to the space ofn x n matrices given by the boundary measurement matrix is Poisson with 
respect to any standard Poisson bracket on &n and the standard Sklyanin bracket on GL„. 
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Remark 2.2. Note that the definition of the R-matrix Rafi in [21] contains a superfluous 
factor 1/2. 

Networks in an annulus that we study in this paper are obtained from the above de- 
scribed networks in a disk by a gluing procedure described in detail in Section 5. These 
networks have one source and one sink, both lying on the same connected component of 
the boundary. The other connected component of the boundary does not carry boundary 
vertices, and hence for our networks if'idS) = //"(G U dS), which impHes Tn = T^. 
Therefore, the standard Poisson bracket on Tn is completely described by Proposition 2.2. 



3. COXETER DOUBLE BrUHAI CELLS 

We start this section with describing a particular instance of the Berenstein-Fomin- 
Zelevinsky parametrization [2, 15] in the case of Coxeter double Bruhat cells in GL„. 

Denote S[p^-^ = SpSp+i ■ ■ ■ Sg-i for I < p < q < n and recall that every Coxeter element 
V e 5„ can be written in the form 

(3.1) v = .V[4_,,y 

for some subset / = {1 = i'o < j'l < ... < ik = £ [1, «]. Besides, define L = {I = Iq < h < 
... < In-k = n} by {/i < . . . < In-k-i} = [1, n] \ /. 

Lemma 3.1. Let v be given by (3.1), then 

Proof. We use induction on n. Denote the right-hand side of the above relation by v. The 
index n - 1 belongs either to / or to L. In the latter case l„-k-i = n - 1, and we have v = 

S\k-,.n] ■ ■ ■ *[ii,i2]«[l,il] = «[in,n-l] ■ ■ ■ ^HhJ2]^Hl.h] '^"-1 and V = ^5„_i.S'[/„_j_,,/„_j_,j ■ ■ ■ Z,].?!!,;,]. 

Then v = v's„-i and v = s„-iv', where v', v' are Coxeter elements in 5„_i corresponding to 
index sets / \ {n} U {n - 1} and L \ {«}, and hence vv = v'v' = 1 by the induction hypothesis. 
Otherwise, if n - 1 belongs to /, we interchange the roles of v and v and use the same 
argument. □ 

Lemma 3.2. The permutation matrix that corresponds to a Coxeter element v is 

k n-k 
7=1 7=1 

Proof. We use the same inductive argument as in the proof of Lemma 3.1. Assuming that 

n - 1 E L, the relation v = v' s„-\ and the induction hypothesis imply v - (ei,-, + . . . -H 
ek-iix-i +ei^i+ ... + e„-i/„_j_2 + e„„)ien + ... + Cn-in-i + em-i + e„-i„) = eui + ... + eij_,„ -H 
ei^i+ ... + e„-u„_t_:, + e„n-\ as claimed. □ 

Let now (m, v) be a pair of Coxeter elements and 

= {1 = /q < ;| < ... < = n), 
= {1 = i'o < i\ < ... < i^- - n), 

^^■^^ L^ = [l = lt<l\<...< /„V_i < - n}, 

L- = {1 = /q < /i < . . . < < = n) 
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be subsets of that correspond to v and u ' in the way just described. For a set of com- 
plex parameters c^, . ..,c~_^;c\, . . .,c'^_^;d\, ...,d„, define matrices D = diag(<ii, . ..,dn), 

IT-I 

7 J 

(3.3) q = <ea.a^u j e [hk% C" = c'ea^Ua, j e [l,r], 



and 



Z+-1 



C| = ^ c^ea,a+u ye [l,n-k% Cj = ^ c^ea+i,a, ye l,n-k ]. 

Lemma 3.3. A generic element X € G"'^ can be written as 

(3.4) X = (1 - Cp-i •••(!- q-_)->D(l - C^r' •••(!- C[)-\ 
anJ its inverse can be factored as 

(3.5) = (1 + C„V)-' •••(! + C\r'D-\l + Cr)-' •■•(! + C,-)-'. 
Proof. It is easy to see that 

(l-CT)-'=£-_,(cI_,)...£I (cl ), 

(l + q)-' = £^^(-c;^^••£^_l(-c^_^), 

(l + Cp-i=£-(-c- )•••£- (-C- ). 



(3.6) 



Then, by (2.2) and (3.1), a generic X G G"'" can be written as in (3.4). Next, the same 
reasoning as in the proof of Lemma 3.1 implies that 



(i_c[)...(i-q\) 



(i + c„V)-'---(i + cf)-i 



and, similarly. 
Therefore, 



(1 - c;.) • • • (1 - CD = (1 + C7)-> •••(! + q--)-i. 



x-i = (1 - c[) • • • (1 - q^.)o-Hi - q--) • • • (i - cp 

= (1 + c„V)-' •••(! + c[)-iD-Hi + C7)-i •••(! + c^)->. 



The network A^^^,, that corresponds to factorization (3 .4) is obtained by the concatenation 
(left to right) of 2n - 1 building blocks (as depicted in Fig. 1) that correspond to elementary 
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matrices 

E-^_,{c- _,),..., EZ (cZ ),D,e;. (4 ),..., El^icl,), 

V-l t--l k+-i ' 

• ■ ■ ' £4(4)' ■ ■ -^^-i^-iX E\{c\) ■ ■ ■ £^_,(4_,). 

This network has 4(n - 1) internal vertices and 5n - 4 horizontal edges. Similarly to the 
case of generic networks discussed in Section 2.4, one can use the gauge group action to 
change the weights of all horizontal edges except for those belonging to block D to 1 . 

One can use the network Nu,v to derive expressions for factorization parameters c?„ c^, 
c~ in terms of matrix entries of X. These formulas are a simple particular case of the 
general formulas by Berenstein-Fomin-Zelevinsky for restoring factorization parameters 
in double Bruhat cells [2, 15]. 

For a matrix A, denote by Arj',. ; V^' its minor formed by rows numbered n < • • • < and 
columns pi < ■■■ < Pk- 



Lemma 3.4. For any i e [1, n], 



for any i e [l,n - 1], 



x\Y.- 



yii,<-ii ■ 



XI' '1... x" ' 



C; = ._ ._ , ( , < I < I , 

1. 



ii.iT-n 



1 j -j 



+1 [1,17-11 



-r 'J ['''Ji 

iX rt.rt+1 [i.rt-11 

;+ i+ iii+-ii 

/I,...,/"!" 

x.l .i 



I 1+1+1^1 1+1 ' 



•I 1141 

j-i ^'j' 



Proof. Since by (3.4), x[J'!] = = di---di, the first formula follows easily. Next, 

note that, for ij < i < iz^^, there is exactly one directed path in Nu,v that joins the jth 
source with the j'Jth sink. The weight of this path equals cr ^ • • • c~4i- . Thus, there is 
a unique collection of vertex-disjoint paths in Nu,v joining sources i[, i~, i with sinks 

By Lindstrom' s Lennma [24, 27] , this imphes that X._ .'_ is equal to the product 

^ J 

of weights of these paths. Clearly, for i~ < i < j't,,, x'!' 'f' = c'Xl' '1' . Also, 

J J"^^ 'q, — ,'j ' — 

X!' '!!''^ ' - cT_diX'l" '':L . Formulae for c" follow from these relations. Formulas for ct 
are derived in a similar way. □ 

Let us now introduce some combinatorial data that will be useful in the following sec- 
tions. 

Let us fix a pair (m, v) of Coxeter elements, and hence, fix the sets given by (3.2). For 
any ; e [l,n] define integers and i^f by setting 



(3.7) ef = 



if i = if for some < j <k± 

1 otherwise 
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and 

(-1 

(3.8) ^r = i{l-sf)-J]el; 

note that by definition, - I, - 0. (Here and in what follows a relation involving 
variables with superscripts ± is a shorthand for two similar relations: the one obtained by 
simultaneously replacing each ± by +, and the other, by -.) Further, put Mf = : a = 
1, . . . , i} and 

(3.9) kf - maxU : i"; < i}. 
Finally, put 

(3.10) £, = £^+£7 

and 

(3.11) Xi = i+l-Y,eii- 

13=1 

Remark 3.1. It is easy to see that there exist distinct pairs (u, v) and (u', v') that produce 
the same n- tuple £. The ambiguity occurs when £, = 1 for some i € [2,n- 1]. By (3.10), 
this situation corresponds either to £t = 1, £^ = 0, or to £t = 0, £^ = 1. Consequently, the 
number of pairs {u, v) with the identical n-tuple e equals 2 power the number of times £, 

takes value 1 . 

Lemma 3.5. (i) The n-tuples s'^ = (sf) and = (^^) uniquely determine each other. 
(ii) For any i € [1, n], 

J* /or some < j <k± 
otherwise 



Z4 

j8=l 



Mf ^[kf-i + i,kf]^[i-J]4J-Yj 4^- 

Proof, (i) Follows form the fact that the transformation £* i-> ^* defined by (3.8) is given 
by a lower-triangular matrix with a non-zero diagonal, 
(ii) By (3.8), the first equality is equivalent to 

(3.12) i^=j + Y,4- 

j8=l 

By (3.7), the latter can be interpreted as counting the first elements of (£*): exactly j of 
them are equal to 0, and all the other are equal to 1 . 

The second equality follows trivially from (3.7) and (3.8). 



(iii) For any j e [1, n], 



j ifi ■■ 



kf = i 



(iv) For any i € [l,n], 
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(iii) For ; = i"^, follows immediately from (3.9) and (3.12). For i + i^., the same counting 
argument used in (3.12) gives 

,8=1 

(iv) Follows from parts (ii) and (iii). □ 

Remark 3.2. (i) If v = s„_i ■■■ s\, then X is a lower Hessenberg matrix, and if m = 
s\--- Sn-\, then X is an upper Hessenberg matrix. 

(ii) If V = Sn-\ ■ ■ ■ s\ and u - s\ ■ ■ ■ Sn-i, then G"''' consists of tridiagonal matrices with 
non-zero off-diagonal entries (Jacobi matrices). In this case = I~ = [l,n], = 1 and 

= for J = 2, . . . , n. 

(iii) If M = V = s„-i ■ ■ ■ si (which leads to - [1, n], = {1, n|), then elements of G"'^ 
have a structure of recursion operators arising in the theory of orthogonal polynomials on 
the unit circle (see, e.g. [32]). 

(iv) The choice u - v - (sisj, ■ ■ •)(-^2-?4 • • •) (the so-called bipartite Coxeter element) 
gives rise to a special kind of pentadiagonal matrices X (called CMV matrices), which 
serve as an alternative version of recursion operators for orthogonal polynomials on the 
unit circle, see [6, 32]. 

Example 3.1. Letn = 5, v = 53^152 andw = i3i2.s1.s4- The network Ar„_v that corresponds 
to factorization (3.4) is shown in Figure 3. 




Figure 3. Network representation for elements in (7*3*2*1 i4,*4i3*i 



A generic element X e G"'" has a form 



d\ 


xiic| 


X12C2 





\ 


c\xn 


d2 + Cj X12 


X22C2 








CIX21 


^2-^22 


ds + C2X23 


dsCj 





C3X31 


C3X32 


C3X33 


c/4 + X34 













C^d4 


ds + C4X45 , 



One finds by a direct observation that fc"^ = 3 and /"^ = {/q, /j, /j} - {1,3,4,5}, 
and hence L"^ = {l^,l'^,l2) - {1^2,5}. Next, m"' = ^451^2^3, therefore, k~ = 2 and 
/" = {«o,«i ' ~ {1.4,5}, and hence - {Zq, /j^, , '3 } = {1,2, 3,5). Further, 

£^ = (1,1,0,0,0), £- = (1,1,1,0,0), 

and hence 

C = (0,-1,1,2,3), r = (0,-1,-2,1,2). 
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Therefore, 

(Kti = (0,0,1,2,3), (A;-)ii = (0,0,0,1,2), 

and hence 

(Mt)U = ([0, 0], [-1, 0], [-1, 1], [-1, 2], [-1, 3]), 

(M;)U = ([0, 0], [-1,0], [-2, 0], [-2, 1], [-2, 2]). 

Finally, e = (2, 2, 1, 0, 0) and x - (0, -1, -1,0, 1). By Remark 3.1, there is one more pair 
of Coxeter elements such that produces the same 5-tuples: v' = S4S1S2S3 andw' = S2*i*3*4- 
To illustrate Lemma 3.4, we find C3 and c^. First, I = < 3 < i'[ = 4, so the 
corresponding formula in Lemma 3.4 gives = X^/Xy Second, 4 = /j", so the other 
formula in Lemma 3.4 gives = 'Ij/Xj^lJ'^j. It is easy to check that the right hand 

sides of both formulas indeed produce correct answers. We wiU use this example as our 
running example in the next section. 

We conclude this section with a proposition that explains how to recognize an element 
of a Coxeter double Bruhat cell in GL„. For any two subsets R,P c [l,n] and a matrix 
X € GL„ we denote by X{R, P) the submatrix of X formed by rows r e R and columns 
peP. 

Proposition 3.1. An element X 6 GL„ belongs to a Coxeter double Bruhat cell if and only 
if the following conditions hold for any / e [1, n - 1]." 

(i+)rankX([/+ l,n],[l,/]) = 1; 

(L) rankX([l, /],[/+ l,n]) = 1; 

(ii+)rankX([l, «],[!,/])> 1 implies X{[1+ 1, «], [1, / - 1]) = 0; 
(ii_)rankX([l,/],[l,«]) > 1 implies X{[\, I - 1],[/+ l,n]) = 0. 

Proof LetZ e G"-". Note that for any p,r e [l,n], the rank of the submatrix X([r,n], [l,p]) 
does not change under right and left multiplication of X by elements of B+. Since G"'^' c 
B+mB+, this means that we only need to check conditions (i+) and (ii+) for the permutation 
matrix u, for which it is clearly true in view of Lemma 3.2. Similarly, conditions (i_) and 
(ii_) reduce to considering v. 

On the other hand, let X satisfy condition (i+) for any I e [1, « - 1] and let i'[ be the 
largest index such that x,- 1 0. Condition (i+) for / = 1 implies i^ > \. Further, condition 
(i+) for ; = - 1 implies X([;Y + 1, n], [1, ~ 1]) - 0. Similarly, we can define to be the 
largest index such that x,-,- and conclude from condition (i+) for / = that i^ > ij, 
and from condition (i+) for I - i^-l that Xdi^ + 1, n], [i], - 1]) = 0. Continuing in this 
manner, we construct a sequence I- = {I = i^ < i^ < ... < = n] such that 

(3.13) Xi-i-_^ ^ 0, X([j; + 1,«], [j~ 1, i: - 1]) = 

for s e [ 1 , fc" - 1 1 . Multiplying X on the right and on the left by appropriate elements of 
B+, we can reduce it to a matrix X' - (x^ ) satisfying (3.13) and such that x'._ - da-, 
x'^j = Moreover, condition (ii+) imply that x'-j = for i^-i < j < i < is. To 

summarize, the lower triangular part of X' is + e,-,- + . . . + e,- 

Now, let s be the smallest index such that i^ > s + 1. Then i^ - 2,...,i^_^ - s. 
Consider the (s + l)st column of X'. Entries xj^^j are zero for i > i^ due to (3.13) and 
for 1 < J < i~ due to the properties of X' described above. Since X' is invertible, this 
means that xi 5+1 0. Then the right multiplication by an invertible upper triangular matrix 
reduces X' to a matrix X" such that x"^_^^ = 1 and the rest of the first row entries are 
equal to zero, while the lower triangular part of X" and the entries in the strictly upper 
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triangular part made by previous reductions are left unchanged. Comparing with Lemma 
3.2, we see that the lower triangular part and the first s rows of X" coincide with those of 
a permutation matrix corresponding to some Coxeter element v of 5„. Continuing in the 
same fashion, we can eventually reduce X" through right multiplication by upper triangular 
matrices to the permutation matrix u, thus showing that X e B+mB+. 

The same argument can be used to show that X e B_vB_ for some Coxeter element v, 
based on conditions (i_) and (ii_). This completes the proof. □ 



4. Inverse problem 

4.1. In this section we show how an element X of a Coxeter double Bruhat cell G"'^ that 
admits factorization (3.4) can be restored from its Weyl function (1.3) up to a conjugation 
by a diagonal matrix. 

Recall various useful representations for the Weyl function m(A; X): 

(4.1) miA;X) = iiAl-Xr'eue,)=^ = A 

Here e, denotes the vector (5,a-)^^i of the standard basis in C", ( •, •) is the standard iimer 
product, piA) is the characteristic polynomial of X, q(A) is the characteristic polynomial of 
the (n- l)x(n-l) submatrix of X formed by deleting the first row and column, and 

hjiX) = iXj)n = (XJeuei), ; e Z, 

is the jth moment of X. (Only moments with nonnegative indices are present in (4. 1), how- 
ever, hj{X) for j < 0, that we will need below, are also well-defined, since X is invertible.) 
In what follows, when it does not lead to a confusion, we occasionally omit the argument 
and write h j instead of hj(X). 

To solve the inverse problem, we generahze the approach of [14], where only the cases 
of symmetric or Hessenberg X were treated. The main idea stems from the classical mo- 
ments problem [1]: one considers the space C[A,A^^]/ det(A - X) equipped with the so- 
called moment functional - a bi-linear functional { , ) on Laurent polynomials in one vari- 
able, uniquely defined by the property 

(4.2) {A',Aj)=hi^j. 

X is then realized as a matrix of the operator of multipUcation by A relative to appropriately 
selected bases {p'l'iX)]"!^, {pjW)'i=o bi-orthogonal with respect to the moment functional: 

{p-iA),p^{A)) = 6ij. 

For example, the classical tridiagonal case corresponds to the orthogonalization of the 
sequence 1 , /I, /I" ' . Elements of G*"-' "*''*"-' ' *' (cf. Remark 3. 2(iu)) result from the 
bi-orthogonaUzation of sequences 1, A, ... , A"~^ and A~^, A^~", while CMV matrices 
(Remark 3.2(iv)) correspond to the bi-orthogonaUzation of sequences I, A, A~^,A^, . . . and 
\,A-\A,A'^,.... 

For any / e Z, «' e N define Hankel matrices 

and Hankel determinants 

(4.3) A® = det-TY,®; 
we assume that Ag = 1 for any / e Z. 
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Remark 4.1. (i) Let X be an n x n matrix, then it follows from the Cay ley-Hamilton 



theorem that, for ; > n, the columns of ' are linearly dependent and so ' = 
(ii) In what follows we will frequently use the identity 



(4.4) 



A®, A®, - 



= A<'-^'Af^'>-(Af) , 



which is a particular case of Jacobi's determinantal identity. In particular, for i = n, (4.4) 
and the first part of the Remark imply A*,'"'^ aJ,'""'^ = (a*,'^) for any /. 

The main result of this Section is 

Theorem 4.1. IfX e G"'^ admits factorization (3.4), then 



di 



(4.5) 



for any / e [l,n]. 



A^^'+^A*"'-'' 
I 1-1 



1-1 



i+l 



i+i 



A^- 



i+l 



Ca' 



(x,_,+i) 
i-i 

^ ^1-1 



Remark 4.2. Formulae (4.5) allow us to restore an element X e G"'" from its Weyl function 
m{A; X) only modulo the diagonal conjugation. Indeed, it is clear from (1.3) that m(A; X) = 
m{A; TXT'^) for any invertible diagonal matrix T = diag(fi, . . . , f„). On the other hand, 
under the action X h-> TXT \ factorization parameters di,cf in (3.4) are transformed as 
follows: di (-> di, cf i-> {tilti+iY^cf, thus leaving the left-hand sides in (4.5) unchanged. 

4.2. The rest of the section is devoted to the proof of Theorem 4.1. The proof relies on 
properties of polynomials of the form 







hi-M 


^/-i+2 • 


■ h+i 


(4.6) 


'pf{A) = det 


hi 


hi+i 


■ k+i 






1 


A 


■ A' 



To prove the first equality in (4.5) we need two auxiliary lemmas. 

Lemma 4.1. Let m e [1, n - 1] and X„ be the mx m submatrix ofX e G"'^ obtained by 
deleting n-m last rows and columns. Then 

(4.7) ha{X„) = haiX) 

for a e [Xm — m + l,x,„ + m]. 

Proof. It is enough to prove the claim for X G G"'" that admits factorization (3.4). It is clear 
that X„ does not depend on parameters c*, . . . , c*_j, dm+i, ...,d„. Moreover, X^ e G"""^", 
where m„ and are obtained from u and v, respectively, by deleting all transpositions s, 
with i > m. Consequently, the network Nu„,v„ can be obtained from the network A^^^v by 
deleting all the edges above the horizontal line joining the mth source with the mth sink. 
Note also that if a > then ha(X) is the sum of path weights over all paths from the first 
source to the first sink in the network obtained by the concatenation of a copies Nu,v Thus, 
hai^m) = haiX) as loug as none of the paths involved reaches above the mth horizontal 
level. The smallest positive power of X such that in the corresponding network there is a 
path joining the first source to the first sink and reaching above the mth horizontal level is 
r = r* + r~, where = Tmn{j : i'J >m + 1}. By (3.9), r* = k^ + l. Therefore, (4.7) holds 
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for a G [Q,k^ + + !]• Lemma 3.5(iii), (3.10) and (3.11), the latter interval coincides 
with \Q,Xm + m\. 

Next, consider the network Nu-^,v-^ that represents X"' corresponding to factorization 

(3.5) . Note that this network differs from Nu-\v-^ ■ In particular, in N all "north-east" edges 
are to the left of any "south-east" edge. Once again, the network ^„^i ^,-\ is obtained from 
the network Nu-Kv-^ by deleting all the edges above the horizontal line joining the mi\\ sink 
with the mth source. The smallest positive power of such that in the corresponding 
network obtained by concatenation of copies of Nu-\v-^ there is a path joining the first 
source to the first sink and reaching above the mth horizontal level is ? = r+ + r - 1, where 

= min{7 \ >m+\} and = minfj : /J > w -H 1}. The difference in the formulas for r 
and r stems from the difference in the structure of the networks A'^ and N: the latter already 
contains paths from the first source to the first sink that reach above the first horizontal 
level. Consequently, it is possible that r = 1 for some m > 1, whereas r > 1 for any m> I. 

One can define combinatorial parameters sf and kf similarly to (3.7) and (3.9) based 
on the sets rather than on 7=^ (cp. (3.2)). It follows immediately from definitions that 
sf = \ - for ; e [2, n - 1] and = = 1. One can prove, similarly to Lemma 3.5(iu), 
that kf = i - e^, which translates to kf = Yj'i3=i ~ ^- Since = fe* -i- 1, we get 
f = Tj'f}=i ~ 1' '^^'^ hence, by (3.11), ? - m - k,„. If r = 1, then Xm - m + I - 0, and the 
interval [Q,Xm + m] coincides with [Xm - m + l,K„ + m\. Otherwise we can concatenate 
up to r - 1 = m - \ - x„ networks iV„-i and hence (4.7) holds additionally for a e 
[Xm-m + \,-\]. □ 

Lemma 4.2. Letm e [l,n - 1], then 

(4.8) det(A-X„) = ^?'^)(a). 

In particular, 

(4.9) di---d„ = 
Proof. Let det(/l - X^) - X" + Y^'^^q Omi'V- Then the Hamilton-Cay ley theorem implies 

m—l 
(=0 

for any a € Z. By Lemma 4.1, this relation remains valid if we replace ha+iiXm) with 
ha+i = ha+i(X) for ; = 0, . . . ,ot, as long as Xm - ni + 1 < a < x,„- This means that, after 
the right multiplication of the matrix used in the definition (4.6) of Pm"^ by the unipotent 
matrix 1 -I- J^^Iq am/se^+i^m+i, one gets a matrix of the form 

"H^^ \ 

lA ■■ A"'-^ det(/l-X„) j' 

and (4.8) follows. Since detX„ = di ■ ■■dm, (4.9) drops out immediately from (4.8) and 

(4.6) after substitution A-Q. □ 

Remark 4.3. Combining Remark 4. 1 (ii) with (4.9) for m = n and taking into account that 
detX - d\---dn, we see that for any / 



a('<„ + 1) 
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which implies that for any / 

(4.10) A® = Ai«-i>detX'-'i-". 

Now, the first formula in (4.5) is an easy consequence of (4.9). To be in a position to 
prove the second formula in (4.5), we first need the following statement. For any i e [1, n] 
define subspaces 

= span{e[, . . . , ej], £j = span{ei, . . . , e,}. 

Besides, put 

!-l 

(4.11) rf = i-lf-'^'^dT'^ Yl "^^'^^ ' ^ [2,n], 

7=1 

where sj" are defined in the proof of Lemma 4.1, and yf = I. 
Lemma 4.3. For any i €[l,n] one has 

(4.12) yleJ=elX^' mod Xjli 
and 

(4.13) yJei = X^'ei mod j. 
In particular, 

£1 = span{e[X^i , . . . , e[X^^ ), £t = span{X^r ei, . . . , X^'ei). 



Proof. A proof for (4.12) was given in [14]. We present it here in order to keep the paper 
self-contained. The case of (4.13) can be treated similarly. 
For any X given by (3.4), consider an upper triangular matrix 

v = D(i- cpr\i - cp_^r' •••(!- c^r'- 

Note that V is the upper triangular factor in the Gauss factorization of X. 
By (3.3), e^Ct = for r < and r > jt, and hence 

e^(l-Cp-' = 



e' mod 



r < 
r > i* 



Thus, for j€[l,k% 

el V = dt*el il-Cpr'---{l-Ctr' 

'j-i j-i 'j-i " 

J-I 'J-l J 'J i 

= flf,+ ct ■■■ ci , el mod £t,. 

.1-1 'j-i 'j 1 'j 'j 1 

A similar argument shows that eJV e £p_i for r < This implies 



(j-i 



■■4 -1 



ej+ mod 



ej+ mod £^*_i. 
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Besides, e\ XV ' - mod £.X , hence the above relation can be re-written as 

'U '1-1 '1-1 



(4.14) e^Xj = 



ej+ modX^+_j. 



On the other hand, for / e - 1], define m > so that / + m -i- 1 is the smallest 

index greater than I that belongs to the index set . Then 

eJV-' = e](\ - CJ) ■ ■ ■ (1 - Ct)D-' mod L],^ 

= {i-ir^'ct ■ ■ ■ c;,„flf,;Li) ^L.i mod 

The latter equality implies 

(4.15) e[y-" = ((-l)'«-^c| • • • cl_^d;^ ■ ■ ■ dj}^ el mod 

for any l'^ e distinct from 1 and n. Note now that = i if and only if a = s/" (this 

can be considered as an analog or (3.12)). Furthermore, dJ^.^ ■ ■ ■ dj} = W'j^^i dj ' . Thus, 
one can re-write (4.15) as 

e{V-^M = f\c^dji^ej mod X^i- 

7=1 

Together with eJVX~^ = ej mod £1^^ this leads to 

= (-1)^-1 ncjrf^f'ef mod £U. 

Combining this relation with (4. 14) and Lemma 3.5(ii), one gets (4. 12). □ 

Example 4.1. We illustrate (4.12) using Example 3.1 and Fig. 3. If 7 > then to find i 
such that e^X^ = jfej mod ^ it is enough to find the highest sink that can be reached 
by a path starting from the source 1 in the network obtained by concatenation of j copies 
of Nu,v Thus, we conclude from Fig. 3, that 

e\X - c/ic^cjgg mod e\X^ = diCiC2d3C^el mod X3, 

e\X^ - d\c\c2d2C^d/[c'^e^ mod X^- 

Similarly, using the network NurKv-^ shown in Fig. 4, one observes that e\X~^ = -c^d^^e^ 
mod £[. These relations are in agreement with (4.12). 

Define Laurent polynomials 

Corollary 4.1. (i) One has 

elpt(X)^eJ, pJ(X)ei^ei, ie[l,n]. 

(ii) For any eigenvalue A ofX, the column-vector (f f (A))"^^ and the row-vector (p~ (A))"^^ 
are, respectively, right and left eigenvectors ofX corresponding to A. 
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Figure 4. Network A^j,-i,v-i for the double Bruhat cell G"'^ from Example 3.1 



Proof, (i) We will only give a proof for pf(A). By Lemma 4.3, or = 1, . . . , j - 1, 

form a basis of hence, taking into account (4. 12), we get 



+ r T 



for some coefficients iXa- By Lemma 3.5(iv), this can be re-written as 



Jte^ 



where either ftj - 1 (if = 0), or tti - 1 (if s"^ - 1). Define a polynomial p(/l) = 
2L=i ^a^"^^ ■ By Lemma 4.3, vectors X°'ei, a g span the subspace X,^;. Therefore, 
by Lenmia 3.5(iv), = for or € [A:" j - i + 2, A:" J. This system of 

Unear equations determines p{X) uniquely as 

(-!)('-')< (t^+t- ,-;+!) 
j-1 



which by (3.11) and Lemma 3.5(iii) gives 



p{X) = 



v('</-l) 



PliX). 



It remains to notice that y^ef = e^X*^' '^'/^(X), and the result follows. 

(ii) Let z'* = (2/)"=! ^ right eigenvector of X corresponding to an eigenvalue A. Then 
elpf(X)z'^ = efz-*, which means that zf = p+(/l)e[z^ = pt(A)zj. Therefore, z^ + 0, 

p\{X) = 4 and (p/'(/l))'Li is a right eigenvector of X. The case of {p^ {Xjf.^ can be treated 
in the same way. □ 

Remark 4.4. The statement of Corollary 4.1i is equivalent to saying that Laurent poly- 
nomials /?f (/l) form a bi-orthonormal family with respect to the moment functional (4.2) 

obtained by the Gram process applied to the sequences 1 , i^i* ,A^i ,.. . and 1 , /l^i , , 

Indeed, 

(p^(^),p-(A)) = e[p;(X)pT(X)ei = (e,-,e,) = d,^. 
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For I G [2, n] define 
(4.16) 



!-l 



i=2 j=l 

where sj = + ej. j 6 [1, nj, similarly to (3.10). 
Corollary 4.2. For any k€Z, 



(4.17) 

Proof. By Lemma 4.3, 



[1,'] 



(=2 'i 'i ) 

By Lemma 3.5(iii), (iv) and (3.11), the determinant in the last expression above is equal, 
up to a sign, to t^!"^''\ By Lemma 3.5(ii), the sign is determined as 



j=2 



7=2 



On the Other hand, (4. 11) implies that nU^rrr = (-l)^'='*'"^'^'r,, and (4.17) follows. □ 

Finally, we can complete the proof of Theorem 4.1. To prove the second relation in 
(4.5), observe that by Lemma 4.3 and Corollary 4.1(i), 

Since by (3.8), (3.10), (3. 1 1) and Lemma 3.5(iii), + A;7 - j + 1 = ki - (i - l)e+, the above 
equality gives 



and so 



Consider the ratio 



Taking into account (4.1 1), we obtain 



A- 



which together with the first relation in (4.5) gives the second one. 
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5. Cluster algebra 

5.1. Let Nu,v be the network associated with X e GL„ and the factorization scheme (2.2). 
We will now construct a network ^ in an annulus as follows: 

(i) For each ; e [\,n\, add an edge that is directed from the i\h sink on the right to the 
ith source on the left in such a way that moving from the ith source to the ith sink in Nu,v 
and then returning to the ith source in along the new edge, one traverses a closed contour 
in the counter-clockwise direction. These n new edges do not intersect and to each of them 
we assign weight 1. 

(ii) Place the resulting network in the interior of an annulus in such a way that the cut 
(as defined in Section 2.4) intersects n new edges, and the inner boundary of the annulus is 
inside the domain bounded by the top horizontal path in Nu,v and the nth new edge. 

(iii) Place one source and one sink on the outer boundary of the annulus, the former 
slightly to the right and the latter sUghtly to the left of the cut. Split the first (the outermost) 
new edge into three similarly directed edges by adding two vertices, a black one slightly 
to the right and a white one slightly to the left of the cut. Add an edge with weight w,„ 
directed from the source to the new black vertex and another edge with weight Wom directed 
from the new white vertex to the sink. 

It is important to note that the gauge group is rich enough to assure the possibility of 
assigning weights as described above, with unit weights at prescribed edges. 

Example 5.1. The network A^°^ that corresponds to A^j,,v discussed in Example 3 . 1 is shown 
in Figure 5. 





'C4 



'C2 



Figure 5 . Network N°^^ for Nu,v from Example 3 . 1 
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Now let X be an element in a Coxeter double Bruhat cell G"''', Nu,v be the network that 
corresponds to the factorization (3.4) and ^, be the corresponding network in an annulus. 
Then A^°^ has 2{n - 1) bounded faces /o,, /i„ i & W,n- 1], which we enumerate as follows: 
each face /o,- contains a piece of the cut and each face /i, does not, and the value of i is 
assigned according to the natural bottom to top order inherited from Nu,v There are also 
three unbounded faces: two of them, adjacent to the outer boundary of the annulus, will 
be denoted /oo,/io> where the first index is determined using the same convention as for 
bounded faces. The third unbounded face is adjacent to the irmer boundary. It will be 
denoted by /o„. 

Recall that faces of A'" ^ correspond to the vertices of the directed dual network (A'" ^,)* 
(as defined in Section 2.4). To describe adjacency properties of iN° yy, let us first consider 
iimer vertices of A^°,,. There are altogether 4n - 2 iimer vertices. For every i e [1, n - 1], 
the fth level contains two black and two white vertices. One of the black vertices is an 
endpoint of an edge directed from the (; + l)th level; it is denoted v^(i). The other one is 
an endpoint of an edge directed from the (i - l)th level (or from the source, for j = 1); it 
is denoted v^ii). Similarly, white vertices are start points of the edges directed towards the 
(i + l)th and the (; - l)th levels (or towards the sink, for / = 1); they are denoted v'^{i) and 
v~{i), respectively. The nth level contains only v,^(«) and v^,(n). 

Now we can describe faces /o,-, / e [0, n|, and /i,-, / e [0, n - 1], by listing their vertices 
in the counterclockwise order Below we use the following convention: if a vertex appears 
in the description of a face with the exponent 0, this means that this vertex does not belong 
to the boundary of the face. With this in mind, we obtain 

fu = (v^(Ov;;(0''^v^(0''' v:,(Ov;(i + i)v:,(/+ i)^-v^(! + i)^-v-(i+ 1)), 

h = {K(i)vt(i)'-K(if'yb(i>w(i + ^Hd + ir'^'v;(i + l)<'v+(i + 1)) 
for ( e [2, « - 2] and 



/lO 

/oo 

/ll 

/oi 

/ln-1 

/o«-i 
fan 



source v^(l)v+(l)v^(l)v,,(l) sink), 

sink v;,(l))v^(l) source) , 

v,(l)v:,(l)v;(2)v:,(2)^Jv,(2)^2v,-(2)), 

v:,( l)v;( l)v-( l)Vfc ( l)v;,(2)v, (2)^^ v:,(2)^J (2)) , 

v,(n - l)v-(n - If'-'vlin - iT'-^vHn - l)v;(n)v;(n)) , 

vlin - Dvlin - lf'-^v-„{n - If'-^in - l)v;(«)v^(«)) , 

vt(n)v„(n)). 



Example5.2. Vertices and faces of A?°^ from Example 5.1 are shown in Figure 6. Consider 
face /i3. As we have seen before in Example 3.1, = = = and Sj - 1, so the 
above description yields fn = (v^(3)v;(3)v^(3)v^(4)v:(4)v-(4)v;(4)). 

In our description of iN"^)* below we use the following convention: whenever we say 
that there are or < edges directed from vertex / to vertex /', it means that there are \a\ 
edges directed from /' to /. It is easy to see that for any i e [0, n - 1], faces /o, and fu have 
two common edges: v~ii + 1) -> v^(0 and v^ii) -> v'^ii + 1). The startpoints of each one 
of the edges are white, the endpoints are black, and in both cases face fu lies to the right of 
the edge. This means that in iN"^)* there are two edges directed from /o, to fu. Similarly, 
the description above shows that iN° yy has 
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(i) 1 - Si+i edges directed from fn+i to fn, 

(ii) 2 - Si+i edges directed from fu to /o,+i, 

(iii) 1 - e,+i edges directed from /o,+i to foi, 

(iv) e,+i edges directed from i to /o, 

for i e - 2], one edge directed from /i„_i to /o„ and one edge directed from /o„ to 
/o„_i. Finally, for \i - 7I > 1 and a,b € {0, 1}, vertices fahfbj in V^l,vy '^''t connected 
by edges. 

Next, we associate with every face /j, in Af° ^ a face weight y,,. We will see below that 
jst are related to parameters cf, di via a monomial transformation. At this point, however, 
let us examine the standard Poisson bracket on Tn",.- As was explained in Section 2.4, 
this bracket is completely described by Proposition 2.2, which together with the above 
description of {N°.^)* impUes the following Poisson relations for face weights: 

{yoi,y\i} = 2yoi)'ii, {y\i,yu+{\ = -(i - SM)yuyu+u 

(5.1) {yii,yoi+\) = (2 - Si+i)yiiyoM, {youyoi+i} = -(1 - Sj+Oj'oiyoi+i, 
(.vo^yii+i) = -Si+iyoiyii+i 

for i e [2, n - 2] and 

(yocyio) = yooyio, {yio,^!) = 23/i()>'()i, 

(5.2) {yw,yn] = -yioyii, lyocyoi) = -yooyoi, 

tVOn.yOn-l} = JOnyOn-l. tVOn.J'ln-l} = -JOnJ'ln-l. 

and the rest of the brackets are zero. 
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Denote by M(A) the boundary measurement for the network j,, and put Hq - WinWout- 
We have the following 

Proposition 5.1. LetX e G"'^ be given by (3.4), then 

MiA) = Ho ((/II + X)-^eu ei) = -Hom{-A; X). 

Proof. Clearly, M(A) is a power series in A~^ with the coefficient of A~'^ equal to (-1)*^' 
times the sum of weights of all paths from the source to the sink that cross the cut exactly 
A; times. Moreover, the leading term of M(^) is HqA-K Denote M(^) = Zk=o(-'^ff^k'i''''^- 
Since the weight of every path has a factor Hq, computing Ht/Ho is equivalent to comput- 
ing the boundary measurement between the first source and the first sink in the planar net- 
work obtained by concatenation of k copies of Nu,v Therefore, Hk/Ho = (X*ei, ei) = hk, 
and 

CO / oo \ 

MiA) = Ho 2(-l)'M"'"' = Ho Yji-^fe^'^i = Ho ((Al + xr'euei) . 

k=0 V k=0 ) 

□ 

Let Hn denote the space of rational functions of the form Q/P, where P is a monic 
polynomial of degree n, 2 is a polynomial of degree at most n- I, P and Q are co-prime 
and P(0) + 0. 

Proposition 5.2. The space of boundary measurements associated with the network N°^^ 
is dense in 

Proof. Let us first prove that any boundary measurement indeed belongs to By Propo- 
sition 5.1, the roots of P are exactly the eigenvalues of —X, hence the degree of P equals 
to n. Next, since M{A) - T,'^^o(-l)''HkA~'^~'^ , the value of M{A) at infinity equals zero, and 
hence the degree of Q is at most n - 1. 

By Proposition 5.1 and (4.1), the coprimality statement is equivalent to saying that Z and 
its submatrix obtained by deleting the first row and column have no common eigenvalues. 
Suppose this is not true, and 1 is a common eigenvalue. Denote X - X - Al. Then 
detX = ^p'^j = 0, and, by the Jacobi detemiinantal identity, 

v[2,n] Y-^l,n-l] _ Hot V _i_ vP,"] vt^'""!] — H 

[2,n] [2,fi-l] [2,fi] 

Thus either xfi'V'' or xl^'"^ ,, is zero. Assume the latter is true (the other case can be 

[2,n] [l,n-l] ^ 

treated similarly). Consider the classical adjoint X of X. Since X is degenerate, X has 
rank one. Since Xu = X\„ = 0, either the first row or the first column of X has all zero 
entries. Assume the latter is true. This means that every (n - 1) X (n - 1) minor based on 
the last n - 1 rows of X equals zero, and so the n x (n - 1) submatrix of X obtained by 
deleting the first column does not have the full rank. Then there is a non-zero vector w 
with wi = such that Xw = 0. Therefore, w is linearly independent with the eigenvector 
ipfiA))"^^ of X constructed in Corollary 4.1(ii), whose first component is equal to 1. We 
conclude that the dimension of the eigenspace of X corresponding to A is greater than one. 
However, due to Lemma 4.3 and invertibility ofX,ei is a cyclic vector for X, which implies 
that all eigenspaces of X are one-dimensional. This completes the proof of coprimality by 
contradiction. The case when the first row of X is zero can be treated similarly. 
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To prove that P(0) + 0, we denote P{,X) = A" + pn-iA."'^ + ■ ■ ■ + po- Then relation 
QiA) = MiA)PiA) yields 

(5.3) Yj(-iypiHk^i = 0, k>0, 

i=0 

with p„ = 1. Relations (5.3) for k € [0,n - 1] provide a system of Unear equations for 
po, ■ ■ ■ ,Pn-\- The determinant of this system equals //[JA^;""''. It is well-known (see, e.g. 
Theorem 8.7.1 in [19]), that the co-primality of P and Q is equivalent to the non-vanishing 
of A*f"^\ So,po = PiO) can be restored uniquely as (-l)"Ai"VAi""^\ which by Remark 4.3 
is equal to (-1)" detX. It remains to recall that detX - d\ ■ ■ ■ d„ Q. 

The density statement follows easily from Theorem 4.1: given M{A), one builds Han- 
kel determinants (4.3) and makes use of formulas (4.5) to restore X, provided Hq and all 
determinants in the denominator do not vanish. □ 



Remark 5.1. (i) Since po 0> equations (5.3) extended to A: = -1, -2, . . . can be used as 

a recursive definition of i/_2, 

(ii) Since m{-A; X) = where p{A) - li"=o(~l)'Pi'^' the characteristic polynomial 
of -X, the Cayley-Hamilton theorem implies that for any k€Z, 

n n 

Yji-^)'PihM = (Yj^-iyPiX'^'euei) = (X'p(-X)euei) = 0. 

(=0 (=0 

Therefore, 

(5.4) Hk = hkHo 
for any k €Z. 

(iii) Denote Q{A) = q„-iA" + \- qo. Similarly to (5.3) one gets 

n 

(5.5) i-iy^'qj = 2 (-'^yPiHi-j-u j e [0, n - 1]. 

i=j+i 

The following proposition is a particular case of Theorem 3 . 1 in [22] . 

Proposition 5.3. The standard Poisson bracket on Tni, induces a Poisson bracket on 'R„. 
This bracket is given by 

(5.6) {M(A), Mill)] = -{A M{A) - /u M(ji)) ^^^^ ~ 

A — fl 

Remark 5.2. Using Proposition 5.1, one can deduce from (5.6) the Poisson brackets for 
the Weyl function m(A) = m{A; X): 

(m(A) - m(Li) \ 
+ m{A)mUi)\ . 

(The derivation of (5.7) from (5.6) can be found in [ 1 3] , Proposition 3 .) Thus a combination 
of Theorem 2.1 and Propositions 5.1 and 5.3 provides a network-based proof of the fact 
that the standard Poisson-Lie structure on GL„ induces the Poisson bracket (5.7) on Weyl 
functions. This fact plays a useful role in the study of a multi-Hamiltonian structure of 
Toda flows. 
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5.2. To compute face weights in terms of factorization parameters cf,di , we introduce 
new notation that makes formulas (4.5) more convenient. First of all, for any / e Z, j € N 
define, similarly to (4.3), Hankel determinants 



(5.8) 



= det(//„.^.;_,_i):,^^i; 



we assume thatzfp = 1 for any / G Z. It follows from (5.4) that Af - i/^A 
Let us fix a pair of Coxeter elements (u, v) and denote e = (£()"=i 

(5.9) xoi = xoiis) = Af^ xu = xuis) = Af'^'\ Ci=c\c:. 

Then formulae (4.5) become 



(5.10) 



di 



^li-^Oi-l 



C; = 



^Oj-^Ij-I xlf 



XOi+1 



Xli-l 
XOi-1 



We now compute the face weights for N°y. 

_ ^-1 _ „l-ei„6i+i-lv2 e,-2 -e.+i 
yOi - - -^Oi-l-^Oi+l -^li-^li-l-^li+l ' 



(5.11) 



yu = 



Cjdi 



'M+i 



l-e, £,+ 1-1 -2 e,- 2-£,-+ 
'■li-l-^li+1 % %-l%+l 



for i G [\,n- 1] and 

1 ifi 

(5.12) yoo = — = , yw = -^ = :>?oiX-^\, 

tin " 



yOn — d„ — X\nXQn-\XQ^X^^_]^. 



We will re-write (5.11) for j = n - 1 in a shghtly difi"erent way. Recall that s„ = 0. Due 
to (4.10), 



(5.13) 

Thus, (5.11) yields 



xl„x-,'„=A'^:'-'^=A^r'\d^iX)'''--", 



xo„=4''"'=4""VetX)''"- 



■n+l 



,, _ „l-e>.-i „-2 „en-i A(n-1) 
y\n-\ - Xi„_2 XQ„_iXQ^_2n^ 



-1 / 1 



detX 



detX 



Define 
(5.14) 



X = X(e) = (Xi)-"i = I Xoi, Xii, . . . , Xo„-l, xin-i,^;,' 



(«-l) 



^(«-2) 



detXl 



andy = y(s) = (yOl'i = (yocjio, • ■ • ,yo«-i,yi«-i)- Theny,- = Ylf^ixJ', where A = (a^Oji 
is an n X n block lower-triangular matrix with 2x2 blocks: 



A = 



Vi 











> 


u 













•^2 


u 


V3 


























•^n-l 


u 


v„ ) 
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with 



(5.15) 



The matrix A is invertible, since det Vi - 1, i e [l,n - I] and det y„ = -1. 



Remark 5.3. Note that the expression for X2n in terms of face weights is independent of s: 

X2n = 0'003'10)"(3'0l3'll)"~^ • • • CVOn-iyin-l)- 

5.3. Given a pair of Coxeter elements (m, v), we want to define a cluster algebra with the 
compatible Poisson bracket given by (5.6). To this end, we use the strategy developed 
in [20]. The first step consists in finding a coordinate system on such that written in 
terms of their logarithms, the Poisson bracket (5.6) becomes constant. Having in mind 
Proposition 2.1, we require this coordinate system to be given by a collection of regular 
functions on Clearly, i > 0, are regular on and hence so are and A^"^ . 

Besides, it was explained in theproof of Lemma 5.2 thatzf*," andzl*,"' do not vanish on'??„, 
hence {A^r\^ and (zll"^' are regular as well. Consequently, by Remark 5.1(i), Hi are 
regular functions on for i < 0, and hence so are Hankel determinants (5.8) for any / e Z, 
/ € N. Therefore, components of x(e) are regular functions on and they are connected by 
an invertible monomial transformation to face weights y(e) that satisfy Poisson relations 
(5.1), (5.2) of the required kind. Therefore we can use x(e) as an initial cluster. Now, 
following [20] , we have to compute the matrix that defines cluster transformations, based 
on the coefficient matrix of the bracket (5.6). 
Define a 2n x 2n matrix 



(5.16) 



B{e) 



U 






y2 

U V3 







with 2x2 block coefficients given by (5.15). Denote by B{s) the (2« - 2) x 2n submatrix 
of B{s) formed by the first In -2 rows and consider the cluster algebra ffig of rank In -2 
with the initial seed S(e) = (x(e), B(e)), so that x„ e [1, 2n - 2], are cluster variables and 

X2n-\, X2n are stable variables. 

Lemma 5.1. Poisson structure (5.6) is compatible with the cluster algebra 

Proof. Let us first revisit standard Poisson structure on fN^^ described by (5.1), (5.2). It 
is easy to see that in terms of the components of the vector y = yis), this bracket can be 
written as {yi,yj] = lOijyiyj, where the matrix Q = Q(e) = (Wy)?"^! is given by 



(5.17) 



Q = 








u 





y2 
-vl 







u 



Therefore, the matrix of coefficients of the Poisson bracket (5.6) written in coordinates x(e) 
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Note that Q defined by (5.17) is invertible. To see that, observe that the block-entries of 
Q satisfy relations VTUVi = U,i € [l,n - 1], and U'^ = 4I2, which impUes that O can be 
factored as 



(5.18) 



















2 "^n-l^ 



^ 




12 I 












if/ ) 



Therefore, detQ = 1, and hence is invertible. 

To find its inverse A^QT^A, observe that if we define 



( 




I2 






12 








•. 

I 







then A = Q.j'^ + y„ ® £„„ with £„„ = (5m(5j„)"^.^i. We then have 

A^Q-i A = JQFJ^ -J{Vn® E„„) + (yj ® £„„) + (yj ® £„„) iV„ ® E„„) = Bis), 

since by (5 . 1 8), the lower-right 2x2 block of equals - i C/ and yj t/y„ = - f/. 

So, B(e) is non-degenerate and skew-symmetric. Thus we can invoke Theorem 1.4 of 
[20]. According to equation (1.5) in the proof of this theorem, compatibility will follow 
from the condition 5(e)0'' - (D 0), where Z) is a (2n - 2) x (2n - 2) diagonal matrix. Since 
B"'(e) = Q^, this condition is obviously satisfied with D = l2„-2. □ 

Our goal is to prove 

Theorem 5.1. (i) The cluster algebra Me does not depend on s. 

(ii) The localization of the complex form ofjl^ with respect to the stable variables X2n-\, 
X2n is isomorphic to the ring of regular functions on 'R„. 

Proof. First, we will compute cluster transformations (2.4) of the initial cluster x(s) in 
directions (Oi) and (li). The transformed variables are denoted xoi and xu, respectively. 
By (5.16), for J e [1, n - 2] the transformations in question are determined by the matrix 



- 1 2 - Si -2 1 - Si+i 
-Si Si - I 2 Si+i - 2 

Therefore, we have to consider the following cases. 



Si+l 
1 - Ei+l 



Case 1: e, = e,+i = 0. Then by (3.11), Ki+i = Xi + l and = J<, - 1, so xq, = 



IS 



transformed into 



XOi = 



Using (4.4) and (5.4), we re-write the numerator as 

1-1 V ' ' '+1 / '+1 V '-1/ 

I i-l I 1+1 \V '-1/ '-1 '-1 / 



1) ^(«,+2) 
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and so 

(5.19) xo,- = Af:;'Uf-^^ - A^A^;^>. 

Similarly, xu = A^^'^^^ is transformed into 



I 

which can be re-written as 

(5.20) Xxi-A^ /l^.^j • 

Case 2: e, = £,+i = 2. This case is similar to Case 1. We have Ki^\ = Ki - I and 
Xi-i = + 1, hence 

(5.21) xoi- -A. zf.^j 

I 

and 

(5.22) xu = ^ , ^ ' ' = ^T'^^"'-'' - 

Case 3: e, = 0, e,+i = 2. We have ;>i,+i = = - 1, and so xq, is transformed into 

The numerator of the above expression can be re- written as 

= Af (zff-2> {Af'^'^f - Af {AfjlAf^l + Af'^'^Af'-'^fj , 

and so 

(5.23) xoi = Af'-^^ [Af'^'^f - Af'^ {^jlAf^l + Af'^'^Af'''^) . 
On the other hand, 



Case 4: s, = 2, £,+1 = 0. This case is similar to Case 3. We have = Xi-i = ki + 1, 
hence 



(5.25) %= '"^c./' ^ 



(z,+2) 
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and 

\ i-i i+i ) 1-1 (+1 V ' / 



Xii 



(5.26) A 

I V ' / ' V '~i '+1 ' ' / 

Case 5: e, = 1, Si+i = 2. We have Xi+i = ki - \, Xi-\ - Xi, so xq; is transformed 
via (5.21) and via (5.24). 

Case 6: e, = 2, e,+ i = 1. We have Xi+i = Xi, Xi-i =;«,■ + 1, so xq; is transformed 
via (5.25) and xu via (5.22). 

Case 7: e, = 0, e,+i = 1. We have Xi+i = Xi, Xi-i = Xi - 1, so xq, is transformed 
via (5.19) and xu via (5.24). 

Case 8: e, = 1, e,+i = 0. We have = Xi + 1, = x,, so xq, is transformed 
via (5.25) and xu via (5.20). 

Case 9: e, = e,+i = 1. We have Xi+\ = Xi-i = Xi, so xq,- is transformed via (5.25) and 
Xu via (5.24). 

Now, let J = n - 1. In this situation transformations of the initial cluster are determined 
by the matrix 

/ En-i - 1 2 - e„_i -2 1 n-Xn-l \ 
\ -e„_i e„_i-l 2 -1 x„-n )' 
Note that x„ does not exceed n - 1, so the last two elements in the first row are always 
nonnegative, and the last two elements in the second row are always negative. By (5.13), 
they contribute to the corresponding relations A^"^ and^^''"~'\ respectively. Therefore, we 
have to consider the following cases. 

Case 10: e„_i - 0. Then x„^i - x„ - \ - x„_2 + 1. and xo„_i is transformed into 



n-l 



Similarly to Case 1, this gives (5.19) for / = n - 1. 
On the other hand, 

= ' 

n-l 

which gives (5.24) for j = n - 1. We thus see that the transformations in this case are 
exactly the same as in Case 7. 

Case 11: e„_i = 1. Then Xn-\ =x„-l = x„-2, and hence 



\ \ J " n-. 

XOn-1 - - 



^(.„-i) 

n-l 

which gives (5.25) for i = n — \. 
Similarly, 

= -J^l ' 

n-l 

which gives (5.24) for i = n-l. We thus see that the transformations in this case are 
exactly the same as in Case 9. 

Case 12: s„_i = 2. Then Xn-\ = x„ - I = x„-2 - 1, and xo„-i transforms exactly as in 
the previous case. 
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On the other hand, 

n-1 

Similarly to Case 2, this gives (5.22) for j = n - 1 . We thus see that the transformations in 
this case are exactly the same as in Case 6. 

Let (m',v') be an arbitrary pair of Coxeter elements, s' be the corresponding n-tuple 
builtby (3.7) and (3.10). 

Lemma 5.2. For any Coxeter elements u', V, the seed S(e') = (x(£'). B{e')) belongs to 
Re- 
proof. First we will show that, in certain cases, mutations of the seed S(e) transform it 
into a seed equivalent to £(£') for an appropriately chosen e' . These situations are listed 
in the table below. In this table, only the entries at which s and e' differ are specified. In 
the first four rows / is assumed to be less than n-1. The second column describes the 
direction of the seed mutation: under the mutation in direction (s, i), the cluster variable 
Xsi is being transformed. It should also be understood that each mutation is followed by 
the permutation of variables with indices (0, and (1, ;) in the new cluster, which results 
in permuting columns and rows 2i - 1 and 2i in the matrix obtained via the corresponding 
matrix mutation. In particular, if x(e') is obtained from x(e) via the cluster transformation 
in direction {s, i), then x(e) is obtained from x(e') via the cluster transformation in direction 
(l-*,0. 



e 


Direction 


s' 


Si = 0, fii+i = 2 


(1,0 




= 1.^.1 = 1 


Si = 2, Si+i = 


(0,0 


s 




Si = l,e,+ i = 


(0, 


s 


= o,«;„ = i 


Si = 2,s,^i = 1 


(0, /) 




-1-^4, =2 


s„-i = 


(0,n- 1) 


e' , = 1 


S„-l - 1 


(1,«-1) 





We will only provide justification for rows one and five of the table. The remaining 
cases can be treated similarly. If e, = 0, e,+i = 2, let e' be defined by ej = 1, = 1 and 

s'j - Ej for /■ t i, i + 1 Then it is easy to check that the matrix mutation in the direction (1,0 
followed by the permutation of rows and columns 2i - 1 and 2i transforms B(s) into B(e'). 
Note also that when e is replaced by s', the corresponding sequence x = transforms 
into a sequence x' that differs from x only in the component x'. - x, - I. This means, that 
x' = x(e') differs from x(e) only in components xo,(e') = zf^"'"'' = xuis) (cf. (5.23) in 
Case 3) and xi,(e') = A^^'^ - xoi(s). Thus, we see that the seed mutation in direction (1, 
of the initial seed of Me transforms it into a seed equivalent to the initial seed of ■ 

Now consider the case e„_i - 0, j = 1. Then p<„_i = x„-l, x'^_^ - x„-\ - 1 and x'^ - 
Xn-1. The fact that the matrix mutation in direction ( 1 , n - 1 ) followed by the permutation 
of rows and columns 2n - 1 and 2n transforms B{e) into B(e') becomes easy to check once 
we recall that, by(3.11),n-l->i;„is always nonnegative. As was shown in Case 1 1 above, 

xu-i(s) = A^^:f = = xon-iis'). Also xo„-i(e) = ^|f_7'' = Zl^,''''^ = Xi„_i(£'), 

which completes the check. 

To complete the proof of the lemma, it suffices to show that, for any s', the seed I.(s') 
is a seed in J^lgco) for e*"^ = (2, 0, . . . , 0). This can be done by induction on ^"^2 ^'i- Indeed, 
if e^_j i= 0, then 2(e') can be obtained via a single mutation from E(e), where s differs 
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from e' only in the (n - l)st component: e„_i - e^ j - 1 (see the last two rows of the 
above table). Otherwise, if ; e [2,n- 2 \ is the largest index such that £■ 0, then, using 
the table again, we see that 2(e') can be obtained via a sequence of mutations from 
where s = (s[, . . . ,s'.^,e'- - 1 , 0, . . . , 0). The intermediate transformations of the n-tuple 
e in this case are (e\,..., s[_^,s[ - 1, 0, . . . , 0, 1, 0), {e\, . . . , e;_j, e; - 1, 0, . . . , 1, 0, 0), . . ., 

(£;,...,e;_j,e;- 1,1,0,. ..,0,0). □ 

The first statement of Theorem 5.1 follows immediately. We can now drop the depen- 
dence on e in the cluster algebra Jig and denote it simply by 



Lemma 5.3. For any j € [l,n - I], k eZ, the function 
(5.27) x(j,k) 
is a cluster variable in 



Ak) max(0,*:+7+l-2«) 



Proof. Consider the cluster x - x(e) that corresponds to e - (2, 0, ... , 0). In this case 



X2j-\ = ^^^ '^ j £ [1>«], and X2j = ^d^, j € [l,n - 1]. The matrix B{e) can be con- 
veniently represented by a planar graph F, whose vertices are represented by nodes of a 
2x n rectangular grid. Vertices in the top row (listed left to right) correspond to cluster 
variables xi, X3, . . . , X2„-i, and vertices in the bottom row correspond to cluster variables 
X2,X4, . . ., X2„. We will label the jth vertex in the sth row by (s, j), s = 0, 1, /' 6 [2, n] 
(s = corresponds to the top row, and j = 1 to the bottom row). In accordance with (5.16), 
F has edges (i, f) is,j -I) for s = 0, 1 and any 7 e [2, n - 1], edges (0, n) ^ (0, n - 1), 
(1,«- 1) ^ (l,n), (l,n- 1) ^ (0,n), (l,n) ^ (0, n) and double edges (0,;) ^ (I,;) for 
e [1,« - 1] and (1,;) ^ (O,; + 1) for j €[l,n- 2]: 



(5.28) 



V 

• -6- 





We marked the vertex that corresponds to the stable variable X2n differently, as it plays a 
special role in what follows. In particular, we will occasionally perturb a two-row structure 
of transformations of the graph F by "moving around" the white vertex. Note also that, in 
view of the definition (5.27), the cluster variable associated with the vertex {s, j) of F is 
x(;;s + J - l)forje [1,«- 1]. 

Denote by Tp the cluster transformation in direction p. Let us consider the result of the 
composition 

T = T2n-2 ° T2n-4 o ■ ■ ■ o T4 o T2 o 7'2n-3 o • • • o 7^3 o . 

An application of T\ transforms x\ = Hq into Xi = -^(Hi +^2^) = H2 and the graph F into 




Next, an application of 7^3 transforms X3 into X3 = A^^^ (here we use (4.4) with i = I = 2) 
and the graph F into 

• ^ • ^ ■• — • ^ • 



A 
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Continuing in the same fashion and using on the jth step relation (4.4) with i - I - j, we 
conclude that an application of T2n-3 o ■ ■ ■ oTj oT\ tothe initial cluster transforms F into 

• -s • -« • ■ — • > • 



• • < • ^ ■■■ — • >- o 



with the variable X2i-i replaced with xii-i = for all e [l,n - 1]. Similarly, the 
subsequent application of Tin-A o ■■■ 0T40T2 transforms F into 



• ^ • -6 • ^ ■ ■ — • < • a- • 




and replaces X2i with X2i = for aHi e [l,n - 2]. Finally, T2n-2 transforms X2n-2 = 
zf^"-;^ into 

n-l 



where we used (4.10). The corresponding transformation of the graph F is 




To summarize, T results in the the transformation x(j,s + j - 1) i-> x(j,i + j + 1) for 
s = 0, 1, 7 e [l,n - 1], and in replacing the initial graph F (see (5.28)) with r(F) (see 
(5.29)). Observe also that the subgraphs of F and ^(F) spanned by black vertices coincide. 

Arguing in exactly the same fashion, we deduce that for r = 1 , . . . , n - 2, an applica- 
tion of results in a cluster T^x) with the corresponding graph r'^(F) such that (i) the 
subgraphs of F and 7'''(F) spanned by black vertices coincide and (ii) the white vertex is 
connected by simple edges to vertices (1, n - r), (1, n - r - 1) so as to form a cycUcaUy 
oriented triangle. In particular, the graph associated with T"^^(x) is 

• ■< • -6 • ^ ■ • ^ • -6 • 

o 

Furthermore, the cluster variable in T''(x) associated with the vertex is,j), s = 0, 1, 7 e 
[l,n - 1]), in 7"'(F) is xij, s + 7 - 1 + 2r). This claim reUes on repeated appUcations of 
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relations 

x( 1 , - 1 )x( 1 , + 1 ) = Xj';'-'"" x(\,kf + x(2,k), 

x(j, k - l)x(j, k+l)^ x'^^'^'-'-^xU, kf + x(j - l,k)x(j +l,k), je [2, n - 2], 

x{n -l,k- l)x(n -l,k+l) = xf;""x(n -l,k)^ + x^f ^''•*^'-"^x(n -2,k- 

which, in turn, follow easily from (5.27), (4.4), (4.10). 

The same pattern of transformations for cluster variables remains valid also for r > n - 1 . 
However, the graph associated with ^""'(x) has a form 



• ■< • < • -1^ ■ • -6 • -6 • 




and the graph associated with T''(x) for r > n - 1 has a form 



• -6 • ^ • ^ ■ — • • • 




o < • < • • ^ ■ — • • 

a. 



where multiplicities Or, br are given by = 2(r - n) + I, br - 2{r -n) + 2. Thus we have 
shown that x(,j, 1 + y - 1) is a cluster variable in ^ for any j € [l,n- 1], and / > 0. 

To recover xij, 1 - 7 - /) for j e [1, n - 1], / > 0, we act in a similar way, starting with 
the cluster corresponding to e = (2, 2, . . . , 2, 0) and repeatedly applying a composition of 
cluster transformations 

T2n-3 o • • ■ o 7^3 o Ti o T2n-2 ° ■ ■ ■ o T4 o T2- 

Thus, xij, k) is a cluster variable in J?l for A: e Z \ [2 - j, j -2]. To complete the proof, it 
suffices to notice that by (3.11), the range of possible values for xj is [1 - 7, j - 1], and thus 
for any k e [2 - j, j - 2] there exists a cluster x(s) given by (5.14) such that x{j, k) is one 
of its variables. Therefore, this variable belongs to J?l by Lermna 5.2. □ 

Since x{\,k) = H/^ for A: < 2n - 2, the following statements readily apparent from 
Lemma 5.3. 

Corollary 5.1. For any k < 2n - 2, Hi is a cluster variable in 

To prove the second statement of Theorem 5.1, we would hke to apply Proposition 2.1 
in the situation when V = "Rn, M is the cluster algebra discussed above and x is given 
by (5.14). Clearly is Zariski open in C^", as the complement to the union of hypersur- 
faces ^^""'^ = and po = 0. Condition (ii) is satisfied by construction, and condition (iii) 
follows from Cases 1-12 discussed above. It remains to check condition (i). 

Observe that the ring of regular functions on is generated by 2(n +1) functions 
Po, . . .,Pn-i, ?o> • • ■ . ?n-i. Pq^, (^n"~'^)~^- Clearly, the last two generators belong to 
Recall that coefficients pi satisfy relations (5.3). These relations for k e [-2, n-3] provide 
a system of linear equations, and by Corollary 5.1, the coefficients of this system belong 
to The determinant of the system is ^^""^^ = ^^""'Vpq' '^^^^ vanish on ^?„. 
Therefore, coefficients po,..., Pn-i belong to Finally, coefficients qo,..., qn-i belong 
to J^* due to relations (5.5) that involve po> • • • . Pn-i and //^ for k <n. □ 
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5.4. The cluster algebra ^ built above is tightly connected to the cluster algebra studied 
in [25, 9]. We denote the latter since to get it from J[ one has to fix the values of both 
stable variables X2„-i and X2„ at 1 . Another cluster algebra, an intermediate between J[ and 
J?l2, is obtained by fixing the value of X2„ at 1 ; it is denoted . 

The exchange matrix of ^2 is obtained from B{e) by deleting the last two rows. If we 
take e = (2, 1, . . . , 1, 0) and rearrange the cluster variables as xqi, ;\:on-i> ^11. • • • > -^in-i, 
the exchange matrix will be given by 



where C is the Cartan matrix for A„_i. This gives precisely the initial cluster considered 
in [25, 9]. Other clusters related to Q-systems (in what follows we call them Q-clusters) 
are obtained from the initial one by using the exchange relation, which is identical to (4.4). 
It is shown in Lemma 1.3 in [9] that Q-clusters correspond bijectively to Motzkin paths, 
that is, integer sequences {mi , m„_i ) such that \mj - mj+i \ < 1 . It follows immediately 
from (3.11) that {x\,. . . ,x„_i) is a Motzkin path starting at 0. It is easy to check that 
this gives a bijection between g-clusters corresponding to Motzkin paths starting at and 
clusters x(s) studied above: the former are truncations X2(e) obtained from x(e) by deleting 
the stable coordinates. Any other Motzkin path is a translate of a Motzkin path starting at 
0. The corresponding Q-clusters are described by the following statement. 

Lemma 5.4. Let xi(e) = ((zlf'',zff '^'^)"ri',zf5,""'*) be a cluster in obtained by the 
truncation ofx(e), and r e Z. 



(i) The r-shift x j (e) = {{Af'^'^\ Af'^'^^^^y^^^ , zfj," ''^'^^^ is a cluster in Mi, and its exchange 



matrix coincides with that ofxi^s). 

(ii) Let X2(e) = (Af'\/Sf''^^^)1~^ be the further truncation of Xi(s), then its r-shift x''^{e) = 
{Af'^'^\ Af'^'^^^^'yi^l is a Q-cluster corresponding to the Motzkin path [xi + r, . . . , Xn-i + r) 
and its exchange matrix coincides with that ofxjis). 

Proof, (i) For r - \ and s = (2, 0, . . . , 0) the proof consists in an apphcation of r2n-3 ° 
• • • o Ti to the graph F shown on (5.28) with the white vertex deleted (see the proof of 
Lemma 5.3, and take into an account that X2n = 1 implies via (4.10) that zfjf"^"^'^^ - A^"~^^ 
for any r € Z). To extend this results to any other value e' it suffices to use the cluster 
transformation taking x(e) to x(e')- The case r > \ follows by induction. The case r = -\ 
is treated similarly to the case r - \ with 7'2«-3 o ■ ■■ oTi replaced by 7'2 ° ■ • ■ ° 7'2n-2, and 
the case r < -1 follows by backward induction. 

(ii) Follows immediately from (i). □ 

Consequently, all cluster variables in all Q-clusters (/?a,m„ in the notation of [9]) form a 
subset of [xt{j,k), j e [l,n - 1], k e Z}, where X2{j,k) are obtained from x{j,k) defined 
in (5.27) by setting both stable variables to 1. The correspondence is given by <-> 

X2{a,ma). 

We conclude this section with a proposition that, in Ught of the above fact, impUes the 
central positivity result (Theorem 9.15) in [9]. 

Proposition 5.4. (i) For any s and any j e [ 1 , n - 1], A: e Z, x{j, k) is a Laurent polynomial 
in x(e) with non-negative integer coefficients. 

(ii) For any e and any j € [l,n - I], k,r € Z, xi{j,k) - x{j,k)\x^=i is a Laurent 
polynomial in x\{e) with non-negative integer coefficients. 
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Proof, (i) Define parameters c,, c/, by (5.10). Pick a pair (m, v) of Coxeter elements that 
correspond to e and consider the element X e G" '' defined by (3.4) with factorization 
parameters di and c7 = c,, c"^ - 1. Then Hi - H[)hi(X). This means that for any 
j € [l,n- 1], k € Z, x(j,k) = i/^X2„'^<*''*"'^'"''^^"'A<*'(X), where by Af\X) we mean 

the determinant defined in (4.3) built from hi(X). By Corollary 4.2, A^*' is the prod- 
uct of a Laurent monomial in variables from x(e) 9with coefficient 1) and the minor 

^^/c-jij^l^'^l If - Xj > 0, then, by Lindstrom's lemma, this minor is equal to the sum 

of products of path weights over all collections of non-intersecting paths leading from the 
j lowest sources to the j lowest sinks in the network obtained by concatenating k - Xj 

copies of the network Nu,y. Thus is a polynomial in factorization parameters 

c,, di with non-negative integer coefficients, and the claim follows, since cluster variables 
and factorization parameters are connected by a monomial transformation with no coef- 
ficients. On the other hand, if k - Xj < 0, then, by a well-known determinantal identity, 

[^'"'"'f^l'^^ = (detZ)*-^' l! ^(di--- dnf-'^i and the previous argu- 

ment applies. 

(ii) By Lemma 5.4(i), x[(s) is indeed a cluster in and its exchange matrix corre- 
sponds to B(e). Define parameters cf\d^p by (5.10) with every Hankel determinant .4® 
replaced by ^['"^"^^ Pick a pair (u, v) of Coxeter elements that correspond to s and con- 
sider the element X e G"'^ defined by (3.4) with factorization parameters J, = df^ and 
cj = cf, c+ = 1. Then 7/,+^ = Hrhi(X) for ; e [0, . . . , 2n - 1]. Recursion (5.3) to- 
gether with Remark 5.1 imply that, in fact, H,+r = Hrhi{X) for all j € Z. Therefore, 
xU, k) = Hlx^'-°-''*^^^~^''^Af~''\X), and the rest of the proof is identical to (i). □ 

Remark 5.4. (i) In fact, we can refine Proposition 5.4(i) and prove Laurent positivity of 

x(j, k) with respect to shifted clusters as well. However, this proof needs additional tools 
in cluster algebra theory, and will be pubhshed elsewhere. 

(ii) If factorization parameters in (3.4) are positive, then the matrix X is totally non- 
negative, and so are matrices X* for k = 1,2, .. . and ./Z*./ with J = diag((-l)')"^j for 
k = -I, -2, . . .. This indicates, in particular, a connection between Q-systems and totally 
nonnegative matrices and their network interpretation. This connection is explored in [10], 
Section 7. 

(iii) The quantization of the cluster algebra considered in this subsection is the subject 
of the forthcoming paper [4]. 

6. CoXETER-TODA LATTICES 

6.1. The goal of this section is to establish a connection between the cluster algebra de- 
fined above and transformations of Coxeter-Toda flows. First, consider the Toda hierarchy 
defined by (1.1). Equations on X induce an evolution of the corresponding Weyl function 
m(A;X), which can be most conveniently described in terms of its Laurent coefficients /;,. 
The following proposition is well known in the case of the usual (tridiagonal) Toda flows. 

Proposition 6.1. If X = X(t) satisfies the Lax equation (1.1), then coefficients hi(X) = 
(X'ei, ei) of the Laurent expansion of the Weyl function m{A; X) evolve according to equa- 
tions 

^MX) = hi^kiX) - hiX)hiiX). 
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Proof. If X satisfies (1.1), then so does X\ By rewriting X* as ;r+(X*) + n-(X'') + ;ro(X^), 
we get 



Now, let (m, v) be a pair of Coxeter elements. Coxeter-Toda flows on G"'''/H are in- 
duced by the restriction of the Toda hierarchy to G"'^. To get a more detailed description 
of Coxeter-Toda flows, we choose parameters c, = etc", <i, that correspond to the factor- 
ization (3.4) of a generic element in G" '' as coordinates on the open dense set in G"'''/H. 
Indeed, c„ di are invariant under conjugation by diagonal matrices (cf. Remark 4.2) and 
are clearly independent as functions on C'^/H. 

Lemma 6.1. The standard Poisson-Lie structure on GL„ induces the following Poisson 
brackets for variables Ci, d^: 

(6.1) {Ci,c,+i} = (£,+1 - l)c,c,+i, {di,dj} = 0, {Ci,di} = -Cidi, {ci,di+i] = Cidi+u 
and the rest of the brackets are zero. 

Proof. In view of Theorem 2.1, it is sufficient to compute Poisson brackets for c,, t/, in- 
duced by Poisson brackets (5.1), (5.2) for face weights of the network Af°y. The first 
equation is an easy consequence of the equality y^i = cj^, i e [l,n - 1], (cf. (5.11)) and 
Poisson relations foryo; described in (5.1), (5.2). 

By (5.11), (5.12), yoiy\i = di/di+i for / e [0, n - 1] (here do = 1). Therefore, 

{logdi/dn.i,logdj/dj+i} = {logyoi7i/, logyo;}'! j), i, j e [0, n - 1], 

which equals the sum of the entries of the 2x2 block of in rows 2i + l,2i+2 and columns 
2j + \,2j + 2. By (5.17), each such block is proportional either to U, or to Vk, or to Vj, 
k€. [1, n - 1], given by (5.15). It is easy to see that the sum of the entries for each of these 
matrices equals zero, and hence {<i,/<i,+i , dj/dj+i] = for all i, j e [0, n - 1]. In particular, 
this holds for / = 0, which can be re-written as {d\,djldj+\] - for all j E [0, n - 1]. 
Taking into account that dj = d\{d2ld\) ■ ■ ■ (di/di-i), we get the second formula in (6.1). 



for i €[l,n- 1], j € [0, n - 1], which equals the sum of the two upper entries of the 2 x 2 
block of Q in rows 2i + 1,2; + 2 and columns 27-1- 1,27-1-2. By (5.17), if such a block 
is nontrivial, it is equal either to U, or to V^, or to -Vj, k e [l,n - 1], given by (5.15). 
Since the sum of the two upper entries equals 2 for U and -1 in the other two cases, we 
conclude that {logc,-, logdj+i/dj] = 26ij - - 6ij-i for ; e [1, n - 1], j £ [0, n - 1]. In 
particular, for 7 = one gets {logc,, log di} - -6a for i e [l,n- 1]. Re- writing dj via di 
and Mas before, one gets {logc,, log i/j} = -6ij + Sij-i,i € [l,n-l],7€ [l,n], which 
is equivalent to the last two equations in (6. 1). □ 

Remark 6.1. We could have also computed brackets (6. 1) by specializing general formulas 
obtained in [26] for Poisson brackets for factorization parameters of an arbitrary double 
Bruhat ceU in a standard semisimple Poisson-Lie group. 




i [X' (Z* - 2;r4Z*) - noiX'')) euei) - ^ ((-X* + 2;r_(X*) + noiX>'))x'euei) 
(r^'^euei) - (;ro(X*)ei, ei) (x'euei) = h^dX) - hiXMX). 



□ 



Similarly, 



{logc,-, log J^+i/Jy} = {logl/ci,logdj/dj+i] = {logyoi.logyo^yij}, 
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Due to their invariance under conjugation by elements of H, Hamiltonians Fk(X) = 
jtrX^ of the Toda flows, when restricted to a Coxeter double Bruhat cell G""", can be 
expressed as functions of c,,d„ which, in turn, serve as Hamiltonians for Coxeter-Toda 
flows on G"'^/H. The easiest way to write down F). as a function of c„ J, explicitly is to 
observe that trX* is equal to the sum of weights of all paths that start and end at the same 
level in the planar network obtained by concatenation of k copies of Nu,v In the case ^ = 1, 
we only need to use Nu,v itself, which leads to the following formula for Fx : define /~ and 
r by (3.2) and denote /" U = {1 = j'l, ...,?„, = n], then 

k-\ ii+i 

(6.2) Fi = Flic d)^di+^ ^ (dj + c^-Wj-i + . . . Cy_i • • • Ci,di,) . 

1=1 j=i,+l 

One can use (6.2), (6.1) to write equations of the first Coxeter-Toda flow generated by Fi 
on G"'"/!! as a system of evolution equations for c„ <i,. 

Example 6.1. (i) For our running Example 3.1, /" U /"^ = {1, 3, 4, 5), so (6.2) becomes 

Fi = d\ + d2 + C\di + di+ C2d2 + C2Cidi + c/4 + c^dj, + c/5 + C4J4. 

(ii) Let V = s„_i • • • si, then U = [1, n] and formula (6.2) reads F\(c, d) = d\ + d2 + 
c\di + ... + dn + Cn-idn-i. If, in addition, u - v then e = (2, 0, . . . , 0) and Fi and (6.1) 
generate Hamiltonian equations 

d 

—di = {di,Fx] = {di,Cidi + c,-ic/,_i} = di{cidi - Ci-xdi-x), 
dt 

d 

—Ci = {Ci, Fx] = {Ci,di + di+x + Ci-xdi-x + Cidi + c,+i<i,+i} = c,(<i,+i - di + Ci-xdi-x - cidi). 
dt 

Then a change of variables r2,-i = di, i e [l,n], and ri, = c,<i„ i e - 1], results in the 
equations of the open Volterra lattice: 

d 

—ji = riiri+x - r,-i), i e [l,2n - 1]; ro = r2n = 0. 
dt 

Another change of variables, a, = c,<i?, fo, = <i, + Ci-xdi-x, leads to equations of motion of 
the Toda lattice that were presented in the introduction. Note that a,, Z?, are, resp., subdi- 
agonal and diagonal matrix entries in a lower Hessenberg representative of an element in 
G"'^/H defined by parameters c,, J,. 

(iii) If M = V = s„-x ■ ■ ■ , then s = {2, 1, . . . , 1, 0}, and Hamiltonian equations generated 
by Fx and (6.1) produce the system 

d d 

—di = diicidi - Ci-xdi-x), — c,- = c,(i/,+i - di + Ci+xdi+x - Cidi). 
dt dt 

After the change of variables c; = Cj<i, this system turns into the relativistic Toda lattice 
d d 

—di = di{Ci - Ci-x), —Ci = Ci{di+x - di + c,+i - c,_i). 
dt dt 

Proposition 6.1 combined with Theorem 4.1 suggests a method to solve Coxeter-Toda 
lattices explicitly, following the strategy that was originally applied in [28] to the usual 
Toda lattice. In order to find a solution with initial conditions c,(0),i/,(0) to the Coxeter- 
Toda equation on G" "/!! generated by the Hamiltonian F^, we first define 

~ ^0 

m\X)=m{A;Xm = Yj^i 

i=0 
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to be the Weyl function of any representative X(0) e G"''' of the element in G"'' /H with 
coordinates c,(0), diiO). Let M{A; t) = 2i^o Hi(t)A^''^ be the solution to a linear system on 
Hn described in terms of Laurent coefficients Hi(t) by 

^mt) = Hi^^it), J = 0,1,..., 
at 

with initial conditions //,(0) = For ; < 0, define Hi{t) via (5.3), where (-1)" 'p; are 
coefficients of the characteristic polynomial of X(0). 

Proposition 6.2. The solution with initial conditions c,(0), to the kth Coxeter-Toda 
equation on G"'^/}! is given by formulas (4.5) with hi = hi{t) = Hi{t)/Ho{t), i e Z. 

Proof. An easy calculation shows that h, - hi(f) - Hi(t)/H()(t), i > 0, give the solution to 
the system presented in Proposition 6. 1 with initial conditions /z,(0) - Thus the function 
m{A, f) - YiT=Q '^KO'^ ' ' evolves in the way prescribed by the A;th Toda flow and therefore 
coincides with m(/l;X(f)), where X{f) is the solution of (1.1) with the initial condition 
X{Q). Since coefficients of the characteristic polynomial are preserved by Toda flows. 
Remark 5.1(ii) implies that for i < we also have hi{t) = hi{X{t)). Finally, since the Moser 
map is invertible on G"'^ /H, we see that the system in Proposition 6.1 is, in fact, equivalent 
to the Ath Toda flow on G"'''/H which completes the proof. □ 

We see that for any pair of Coxeter elements (m, v), the Coxeter-Toda flows are equiva- 
lent to the same evolution of Weyl functions. We want to exploit this fact to construct, for 
any two pairs (m, v) and (m', v') of Coxeter elements, a transformation between G"'^/H and 
Qu'y ^jj jjjgj. Poisson and maps the fcth Coxeter-Toda flow into the ^th Coxeter-Toda 
flow. We call such a transformation a generalized Bdcklund-Darboux transformation. The 
term "Backlund transformation" has been used broadly over the years for any transforma- 
tion that maps solutions of one nonlinear equation into solutions of another. To justify the 
use of Darboux's name, we recall that traditionally a Backlund-Darboux transformation 
consists in interchanging factors in some natural factorization of the Lax operator associ- 
ated with a given integrable system. In the case of Coxeter-Toda flows, the same number 
and type of elementary factors appears in the Lax matrix associated with any Coxeter dou- 
ble Bruhat cell. Hence we use the term "generalized Backlund-Darboux transformation" 
even though in our case, re-arrangement of factors is accompanied by a transformation of 
factorization parameters. 

Let us fix two pairs, (m, v) and («', v'), of Coxeter elements and let e = (£()"=!> ^' = 
(£■)"=! be the corresponding n-tuples defined by (3.7), (3.10). We construct a map a^^^ : 
C'^/H — > G" '^' /H using the following procedure. Consider the cluster algebra defined 
in Section 5. Fix a seed E(e) - (x(e),B(e)) in Jl, where x(e) is given by (5.14) and 
Bis) by (5.16). Let be the sequence of cluster transformations defined in the proof of 
Lemma 5.2 that transforms 2(£) into the seed 2(e')- Next, for an element in G"'''/H with 
coordinates c,,d„ consider its representative X g G"'", the corresponding Weyl function 
m(/l; X) and the sequence of moments hiiX), i e Z. Apply transformation t„,v by assigning 
values to cluster variables in the cluster x(e) according to formulas (5.8), (5.9), (5.14) with 
Hi replaced by hiiX). Then apply transformation T§' to x(e) to obtain the cluster x(e')- 
Finally, apply transformation pu'y by using equations (5.10) with e replaced by e' and 
components of x(e) replaced by those of x(e') to compute parameters c'., d'- that serve as 
coordinates of an element in G" " /H. This concludes the construction of cr^^y . 
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Theorem 6.1. The map cr",'^, : G"''7H — > G" /H is a birational transformation that 
preserves the Weyl function, maps Coxeter-Toda flows on G"'^/H into matching Coxeter- 
Toda flows on G" '"'/!! and is Poisson with respect to Poisson structures on G"'^/H and 
Qu'y in induced by the standard Poisson-Lie bracket on GL„. 

Proof. Moments hjiX) are polynomial functions of c;, di for i > and rational functions of 
c„ di for i < 0. Values we assign to cluster variables in x(e) are thus rational functions of 
c,, di. This, combined with the rationality of and equations (5.10), shows that the map 
crjjj,' is rational. It is easy to see that its inverse is cr^J^, which implies birationality. The 

claim that cr" preserves the Weyl function is simply a re-statement of Lemma 5.2, which 
implies that if clusters x(£) and x(e') are obtained from a function M(/l) e according 
to (5.8), (5.9), (5.14), then r|' transforms x(e) into x(e')- The rest of the statement of the 
theorem is a consequence of the invariance of the Weyl function, since Poisson structures 
on G"''7H and G" '' /H induce the same Poisson bracket on compatible with J\ and, by 
Proposition 6.1, Coxeter-Toda flows generated by Hamiltonians on G"'"/!! and G" /H 
induce the same evolution of the Weyl function. 

□ 

To illustrate Theorem 6.1, in the table below we Ust elementary generalized Backlund- 
Darboux transformations that correspond to cluster transformations from a fixed cluster 
x(e) into an adjacent cluster x(e'). The table can be viewed in parallel with the table in 
the proof of Lemma 5.2. Expressions for transformed variables c'j, d'j are obtained by 
combining formulas for cluster transformations with equations (5.10). Variables that are 
not listed are left unchanged. 



e 




Transformation 


Inverse 


£; = 

e,+i = 2 


^; = i 


^' _ Cidi 

d'i = lBk,dl,=di,^+Cidi 


^ «+i 

di = d\(\+c'^, di^X = ^ 


Ei = 2 

£i+i = 




< = dAutf^ ^,+1 =c,+i(l +cO 
dl, = ^,d[=di{\+Ci) 


di,,=d\,,^d\c\,di = ^^ 


ei=l 
e,+i = 


£; = o 


d[=di{\+Ci), = ^ 


r _ « r , _ ^I'^^H 

di+i - d',, + c'd', di — f 

i+l 11^ * d.^^+c.d- 


Si = 2 
Si+l = 1 




< = J,(l+c,), 4,1 = ^ 


c'.d'. c'.^d'.^^ 

di^i = dl, + c'.d'., di = 


Sn-l = 




f,i _ c„_id„_i 
^n-\ ~ d„ 

d'„=dn+ Cn-\dn-l 

J, _ d„d„-i 
n-l d„+c„-id„-i 


^'n-ld'n 

dn-r=d'„_,(l+c'„_,) 

J d'„ 






„, _ c„_id„-i 
'-n-l ~ dn 


^« i4 
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e„_i = 1 



e' , —1 

n-l 



c — 

n-2 dn+Cn-idn-\ 



d'„=d„ + c„-id„-i 



d' 



dndn- 



n-l d„+c„-id„- 



cn-2 = c;_2(i + c;_j) 
4-i = Ci(i + Ci) 

d: 



d„ 



Elementary generalized Backlund-Darboux transformations can be conveniently inter- 
preted in terms of equivalent transformations of perfect networks introduced in [29]. The 
three types of equivalent transformations are shown in Figure 7. Instead of trying to de- 
scribe the general case, we will provide an example. 





Tjfpe 1 





Type 2 



1 



X y 
x' y'f 



1 



W2 



1^ 

1+ Wl M'2 



W2 
l+Wi W2 



y 


X 1 


'y' 


x' 



1+ Wl W2 



1 l^*^'^^ 



Tjfpe 3 



Figure 7. Equivalent transformations of perfect networks 



Example 6.2. Consider the network from Example 5.1. Recall that s = (2, 2, 1, 0, 0) and 

set / = 2. So, £2=2 and £3 = 1, which corresponds to the fourth row of the above table. 
The corresponding transformation consists of the following steps: 

(i) Type 2 transformation with x = vj(3), y = v^(3) and wi = W4 = 1, W2 = , W3 =c^. 

(ii) Type 3 transformation with x = vj(3), y = v~(3), x' = v^(2), y' = v~{2) and wi = C2 , 

W2 = Cj, W3 = ^3, W4 = 1. 

(iii) Type 1 transformation with x = v;J;(2), y = v~i2) and w\ = 02/(1 + C2), W2 = 1, 

Wi - c/2, W4 - 

(iv) The gauge group action at v^(3) that takes the triple of weights (J3, /(I -1-C2), 1/(1 + 

C2)) to (1,^/34/(1 +C2),^;3/(1+C2)). 

(v) The gauge group action at v~{2) that takes the triple of weights (1 + C2, l,Cj) to 
(1,1 + C2,cj;(l +C2)). 

(vi) The gauge group action at v;J;(2) that takes the triple of weights (1 + C2, <i3cj/(l + 

ci), di) to (d2(l + C2), dicyWl + C2)], 1). 

Thus, at the end we have (cj)' = 03/(1 + C2), (cj)' = diC2/[d2il + C2)], and hence 
c'2 = d3C2/[d2il + C2)]. Besides, (cp' = c]" (1 -1- C2), Cj )' = c|, and hence c[ = ci(l -1- C2). 
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1+C2 



Figure 8. Elementary generalized Backlund-Darboux transformation: 
steps (i) and (ii) 



Finally, d'2 - d2(l + C2) and c/3 = d^/il + C2). All these expressions coincide with those 
given in the fourth row of the table. 

Transformations of the relevant part of the network during the first two steps are shown 
in Figure 8. 

Transformations of the relevant part of the network during the remaining four steps are 
shown in Figure 9. 

We can make transformations crj^j,^ more expUcit by using Corollary 4.2. Below we 
write A[i] for the determinant of the leading principal i X i submatrix of a matrix A. Pick an 
element X = X{c, d) € G"'^ defined by (3.4) with factorization parameters J, and c~ = Ci, 

ct = L 

Proposition 6.3. The map crl y transforms coordinates Ci, di on G"-'" jH into coordinates 
c'j,d'. on G" /H given by formulas 

J/ ^ V hi] V Ai-i] 



c'i = Cid] 



[i+l] 







2- 








/ 


V ) 





where s and x (resp. s' and x') are n-tuples (3.10), (3.11) associated with (u,v) (resp. 
(«', v')) and 6k j - k'j - xjfor j e [1, n]. 

Proof. The claim follows from formulas (4.5), (4.16), (4.17) and an easy computation that 
shows that 

r,-ir,+i _ ^2 (^1 • • • ^i-i)^' _ ,2 



(Ji---j,+i)''« 



(-^[i+i]) 



48 



MICHAEL GEKHTMAN, MICHAEL SHAPIRO, AND ALEK VAINSHTEIN 




Figure 9. Elementary generalized Backlund-Darboux transformation: 
steps (iii)-(vi) 



□ 

Example 6.3. Let v = = u' = v' - Sn-i ■ ■ ■ si. Then G"'^/H is the set of Jacobi matrices 
(1.2), which serves as the phase space for the finite nonperiodic Toda lattice, and G"''^7H 
can be viewed as a phase space for the relativistic Toda lattice. Combining Theorem 6. 1 
with Example 6.1, we obtain the following corollary of Proposition 6.3. 

Corollary 6.1. If the entries a,, bi of the Jacobi matrix L evolve according to the equations 
of the Toda lattice, then functions 

_ V /['] \ A'-i] _ V /[t+i] V /[i-i] 

- (Li-')[,(L3-0,_V - {V-%{L^-% 

solve the relativistic Toda lattice. 

Proof. First, observe that an element X featured in Proposition 6.3 is a tridiagonal ma- 
trix whose nonzero off'-diagonal entries are X„+i - di,Xi+u = c,d;. The matrix L as- 
sociated with the same parameters Ci,<i, is related to X via L = DXD~^ where D = 
diag(l,c/i, . . .,di ■ ■ -dn-i). This means that {L!^\i] - {X'^\i\ for any i,k and a, = = 
Cid^. Furthermore, s,- = for i e [2, . . . ,«], Xi - i- l and pcj = for ; e [1, . . . The 
claim then follows from Example 6. l(iii) and formulas of Proposition 6.3. □ 

6.2. It is natural to ask if the classical Darboux transformation X = X-XoX+ D{X) = 
XqX+X- can also be interpreted in terms of the cluster algebra The transformation 
D constitutes a step in the LU -algorithm for computing eigenvalues of a matrix X. A 
connection of the Lf/-algorithm (as well as similar numerical algorithms, such as QR and 
Cholesky algorithms) to integrable systems of Toda type is well-documented, see, e.g. 
[8, 34]. For an arbitrary semisimple Lie group, a restriction of such a transformation to 
a Coxeter double Bruhat cell of type G" " was studied, under the name of factorization 
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dynamics, in [23]. We collect some relevant simple facts about the transformation D in the 
proposition below. 

Proposition 6.4. Let X e N B+. Then 

(i) for any i e Z, hi(D(X)) = hM(X)/hi(X}; 

(ii) /or any u,v€S„, ifX e G"-" then D{X) € G"•^• 

(iii) D descends to a rational Poisson map D : G" ' /H — > G"'^'/H that coincides with a 
time-one map of the Hamiltonian flow generated by the Hamiltonian F(X) = ^ trlog^ X. 

Proof (i) For / > 0, we have hi+i(X) = [X-iXoX^X^yXoX+euei] = di ei) = 

hi{X)hi{D{X)). The case ; < can be treated similarly. 

(ii) It suffices to observe that if Y\ e N- and Y2 e B+ than both statements Y\ Y2 e G"'" 
and Y2Y1 e G" " are equivalent to Yi e B+vB+, Y2 e B uB . 

(iii) Claim (ii) implies that D descends to a rational map from G"''7H to G" "/!!. The 
rest of the claim is an immediate corollary of general results in Section 7.1 in [23]. □ 

For a pair of Coxeter elements (u, v), Proposition 6.4(i) allows us to completely describe 
the action of D on C'^/H in terms of a simple map on !??„. Namely, define rj : Jin 
K by J](MiA)) = AMiA) - Hq. Equivalently, 77 can be described by r] (Y^Zo f^'^''"^) = 
XZo Hi+iA-'-^ . Then Proposition 6.4(i) implies that on G" 7H 

D = pu,v ° x„,v 0770 m„,v, 
where maps p„,v. x„_v. inu,v were defined in the Introduction. 

Remark 6.2. As we have seen in Section 5.4, the shift //, i-> plays an important role 
in the study of Q-systems in [9]. 

To tie together the cluster algebra J?l and the Darboux transformation £), we have to 
descend to the cluster algebra introduced in Section 5.4. we will only need to fix the 
stable variable X2„ to be equal to 1. In view of (5.14), this means that we are dealing 
with double Bruhat cells in S L„ rather than in GL„. In order to emphasize a similarity 
between the classical Darboux transformation D and the generalized Backlund-Darboux 
transformation cj^^^ , we express the former similarly to (1.5). 

Proposition 6.5. D = pu,v ° T'o ° Tu,v, where Tj) is a sequence of cluster transformations 
in 

Proof. Note that in the graphical representation of the matrix B(s) that we employed in 
the proof of Lemma 5.3, passing to the cluster algebra amounts to erasing the white 
vertex and all corresponding edges in the graph F. Consider the cluster corresponding to 
s - (2, 0, ... , 0). By Lemma 5.4(i), the shift Hi Hj+i is achieved by an application of 
T2n-3 o • ■ • o . This means that for v = = s„_i • • • si we can choose r2n-3 o • • • o for 
Tj). Then, for arbitrary pair of Coxeter elements (m, v), can be defined as 

with w = s„-i ■■ ■ si. □ 
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